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OPERATING  CHARACTERISTICS  OF  LOG-NORMALIZER 


FOR  WEIBULL  AND  LOG-NORMAL  INPUTS 


INTRODUCTION 


The  detection  of  the  presence  of  a  weak  signal  of  unknown  location  and 
strength  In  background  noise  of  unknown  strength  is  often  accomplished  by 
comparing  a  candidate  signal-bearing  detection  sample  of  the  observed 
process  with  a  local  estimate  of  the  background  level  based  on  N  samples  of 
the  (hopefully)  noise-only  process.  The  local  neighborhood  can  be  time, 
space,  or  frequency,  depending  on  the  application.  In  order  to  obtain  a 
stable  estimate  of  the  background  level,  the  number  of  samples,  N,  should  be 
large;  however,  if  the  background  is  nonstationary,  nonhomogeneous ,  or 
nonwhite,  or  if  decision  and  processing  time  is  at  a  premium,  N  should  be 
kept  as  small  as  reasonably  possible.  The  tradeoff  between  these  conflicting 
requirements  and  the  dependence  on  the  number  of  normalizer  samples,  N,  is 
of  interest  in  this  study.  Related  work  is  available  in  (1,2,3). 

Since  the  performance  of  the  normalizer  procedure  outlined  above  is 
adversely  affected  by  the  presence  of  any  outliers  or  noise  bursts 
anywhere  in  the  total  of  N+l  samples  used  to  make  a  decision  about  signal 
presence  or  absence  in  the  candidate  sample,  some  form  of  limiting  device 
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should  precede  the  normalization.  The  particular  combination  that  we 
consider  in  detail  here  is  depicted  in  figure  1,  where  in  is  the  natural 


Figure  1.  Log-Normalizer 


logarithm.  The  logarithmic  device  tends  to  saturate  at  large  input 
amplitudes  and  suppress  their  effect  on  the  normalizer  output  z.  The  input 
sequence  of  random  variables,  {x^  ,  (the  detector  output  sequence),  is 
presumed  to  be  statistically  independent  and  limited  to  positive  values, 
giving  logarithmic  output 

yn  *in(xn)  .  (1) 

The  particular  normalizer  we  consider  here  is  described  as  follows: 
call  yQ  the  candidate  signal-bearing  sample  at  the  normalizer  input,  and 

let  y.| ,  y? . yN  be  the  N  noise-only  basis  samples  employed  to 

extract  an  estimate  of  the  background  level  at  the  normalizer  input. 

Despite  the  notation,  these  N  samples  can  (and  probably  will)  surround  the 
candidate  sample  yQ  in  location,  whether  that  be  time,  space,  frequency, 
etc.  The  sample  mean  of  the  normalizer  input  noise-only  samples  is 

N 

*<»>  ■JIi,.  <2> 

n=l 
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while  the  corresponding  sample  standard  deviation  is  defined  as 


<*>  Wn  -  5<v)l^|1/2  • 


(3) 


The  output  of  the  normalizer  in  figure  1,  that  we  consider  here,  is  a 
deflection  measure,  namely 


z  = 


y0  -  v(y) 


3(y) 


(4) 


The  numerator  of  (4)  is  an  estimate  of  the  difference  in  means  (at  the 
normalizer  input)  of  the  candidate  signal-plus-noise  sample,  relative  to  the 
noise-only  samples;  the  denominator  of  (4)  is  a  measure  of  the  inherent 
fluctuation  of  the  background  noise.  The  dimensionless  ratio  in  (4) 
eliminates  the  dependence  on  absolute  levels  in  favor  of  relative  levels. 


The  normalizer  output  z  is  compared  with  a  threshold  T,  and  a  decision 


made  about  signal  presence  in  sample  yQ  according  to  the  rule 


z  >  T:  declare  signal  present  in  y 


*\ 


z  <  T:  declare  signal  absent  in  y 


(5) 


It  is  desired  to  evaluate  the  false  alarm  probability  and  the  detection 
probability,  that  is, 


3 


<: 

v 

3 


4 


3 


n 


Pp  *  Prob  (z  >  T  |  signal  absent  in  yQ), 


PQ  =  Prob  (z  >  T  J  signal  present  in  yQ). 

Both  of  these  probabilities  in  (6)  are  exceedance  distribution  functions, 
that  is,  probabilities  that  random  variable  z  is  greater  than  a  threshold 
value  T.  We  will  be  interested  in  plots  of  (6)  versus  T,  for  various 
signal-to-noise  ratios,  as  well  as  in  plots  of  PQ  versus  P^,  the  latter 
known  as  the  receiver  operating  characteristics. 

The  normalizer  input  random  variables  {y^  are  statistically 
independent  and  identically  distributed,  since  inputs  ![*  have  been 
presumed  to  have  these  properties.  When  signal  is  absent  in  candidate 
sample  yQ,  its  probability  density  function  will  be  taken  identical  to 
that  of  {yn\  1 ;  however,  when  signal  is  present  in  yQ,  its  probability 
density  function  can  be  arbitrary. 


CLASSES  OF  INPUT  VARIABLES 


{*1  N 

xnl  i  t0  lo9~normali zer  in  figure  1 
will  be  taken  from  the  class  of  random  variables  that  can  be  generated  from 
fundamental  independent  identically-distributed  random  variables  {wn] 
according  to  the  rule 

x„  =  a  wj  for  1  <  n  <  N;  a  >  0,  t>  >  0,  w„  >  0  .  (7) 

The  probability  density  function  p^  of  {wn}  is  arbitrary;  the  total  class 
of  random  variables  defined  by  (7)  is  that  yielded  by  allowing  parameters  a 
and  b  to  be  any  positive  constants  (independent  of  n). 

To  fix  this  concept  of  a  class  of  random  variables,  consider  the  case 
where  wn  is  a  random  variable  with  the  fundamental  exponential  probability 
density  function 

Pw(u)  =  exp(-u)  for  u  >  0  .  (8) 

Then  the  exceedance  distribution  Function  of  w  is 

n 

oO 

Qw(u)  =  Prob(w  >  u)  =  dv  Pw(v)  =  exp(-u)  for  u  >  0  .  (9) 

u 


5 
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It  then  follows  from  (7)  that  the  exceedance  distribution  function  of  xn  is 
Qx( u )  *  Prob(x  >  u)  *  Prob(a  wb  >  u)  ■ 


Prob 


(“  >  <u/a),/bj  -  Qw^)'/1>)  -  «p[-(^)'/IJ] 


for  u  >  0 


(10) 


But  this  is  just  the  exceedance  distribution  function  of  a  Weibull  variate* 
with  shape  factor  1/b  and  scaling  (1/a)1^.  Thus,  the  class  of  random 
variables  that  can  be  generated  via  (7)  with  arbitrary  a.b,  from  the 
fundamental  exponential  probability  density  function  in  (8),  is  the  general 
class  of  Weibull  variates,  as  given  by  (10).  (If  b  =  1/2,  x  is  a  Rayleigh 
variate,  for  example.) 


As  a  second  case,  let  wn  be  a  random  variable  with  the  fundamental 
log-normal  exceedance  distribution  function 


where 


0  (u)  =$(-j^n  u)  for  u  >  0  , 

w  * 


(ID 


t  t 

$(t)  s  j*  ds  (2*)'1/2  exp(-s2/2)  s  £  ds  0( s)  (12) 

-*>  -06 

is  the  cumulative  distribution  function  of  a  normalized  Gaussian  random 
variable.  Then  by  an  analogous  procedure  to  (10),  the  exceedance 
distribution  function  of  the  random  variable  xn  generated  according  to  (7) 
is _ _ 

S*e  appendix  A. 


6 
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y  =  >fn(x  )=J?na  +  bjlnw  =  J!n  a  +  b  "  for  1  <  n  <  N  ,  (15) 

n^  n  n  n  ~- 

where  we  define 

v  =  Jfn  w  .  (16) 

n  n 

Let  the  mean  and  standard  deviation  of  be  denoted  by  w(v)  and  o(v), 
respectively.  Then  form  the  normalized  random  variable 


1 
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V  -  w(v) 

v  =  — -  for  1  <  n  <  N  ,  (17) 

n  .  m » 

0(V) 

which  has  mean  0  and  standard  deviation  1.  Substitution  of  (17)  in  (15) 
then  yields  log  output 

y_  *  a  +  B  v  for  1  <  n  <  N  ,  (18) 

n  n  ~  — 

where  constants 

a  »  jn  a  +  b  >j(v)  ,  B  «  b  o(v)  .  (19) 

A  direct  useful  interpretation  of  these  two  constants  in  (19)  follows 
directly  from  (18);  namely,  since  {vj;  are  normalized  random  variables, 

a  =  u(y),  B  =  o(y)  .  (20) 

These  are  fundamental  statistics  of  the  input  to  the  normal izer  in  figure  1. 

Equations  (18)  and  (19)  demonstrate  that  the  output  of  the  logarithmic 

f  1  N 

device  in  figure  1,  for  general  parameters  a,b  and  random  variables  tWpj  1 
in  transformation  (7),  can  be  handled  through  the  linear  transformat  ion  (18) 
of  a  normalized  random  variable,  v  ,  with  zero  mean  and  unit  standard 
deviation.  The  new  general  parameters  a,B  are  given  by  (19)  or  (20),  where 
the  required  statistics  are  mean 
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r  N 

in  terms  of  the  normalized  random  variables  tvn* 1 


As  noted  in  the  paragraph  following  (6),  the  probability  density  function 


of  random  variable  yQ  is  arbitrary  for  the  signal-present  hypothesis. 


Without  loss  of  generality,  let 


yQ  =  n  +  v  v0  . 


where  normalized  random  variable  v  has 

o 


u(vq)  =  0,  a(vQ)  =  1 


The  constants  n  and  v  absorb  the  absolute  scale  of  yQ;  in  fact  (in  analogy 


with  (20)), 


n  =  w(yQ) ,  v  *  o(yQ)  ■ 


When  we  now  combine  (23)  and  (25)  in  the  normalizer  output  z,  as  given 


by  (4) ,  there  follows 


d  rv  -  y(v) 
o 


where  constants 


a 


v 

r  “  0 


WT  '-iMPUWBPVTumirgTa1 
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Thu%  the  general  output  z  of  the  log-normal izer  in  figure  1,  for  the  general 
class  of  inputs  (7),  can  be  expressed  In  the  form  (28)  involving  two 
fundamental  constants  d,  r  in  (29);  an  arbitrary  normalized  random  variable 
vQ;  and  the  sample  mean  and  standard  deviation  of  the  normalized  random 

.  •  |l| 

variables  ^vj 1  according  to  (24). 

A  useful  physical  interpretation  of  the  constants  in  (29)  is  afforded  by 
utilizing  (20)  and  (27),  namely 


p(yQ)  ~  v(y)  o(yQ) 

o(yj  ’  r  =  “STyT  * 


(30) 


Thus,  parameter  d  measures  the  deflection  criterion  at  the  normalizer  input, 

relative  to  the  standard  deviation  for  signal  absent.  The  parameter  r  is  a 

scaling  quantity  reflecting  the  relative  fluctuating  strengths  at  the 

normalizer  input.  The  fundamental  analysis  problem  is  now  to  evaluate  the 

false  alarm  and  detection  probabilities  specified  by  (6),  for  the  output 

random  variable  given  by  (28),  where  d  and  r  are  arbitrary  constants,  and 
r  i  N 

vQ  and  (vnj 1  are  normalized  random  variables. 


i 


1 1 
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CONSTANT  FALSE  ALARM  RATE  PROPERTY 


The  general  output  of  the  log-normal izer  is  given  by  (28).  However,  for 

hypothesis  Hq  where  signal  is  absent  in  candidate  signal -bearing  sample 

yQ,  the  statistics  of  normalizer  input  yQ  are  identical  to  those  of 
f  1  N 

{y  jj,  as  noted  in  the  paragraph  under  (6).  In  this  case,  (30) 
obviously  reduces  to 


d  =  0,  r  =  1  under  Hfl  , 


and  (28)  yields 


v0  -  S(v) 


under  HQ  , 


in  terms  of  the  independent  identically-distributed  normalized  random 


variables  v  and  fv  \ . 

o  l  n'l 


Since  v  in  (17)  is  the  normalized  random  variable  corresponding  to 
n 

logarithmic  distortion  (16)  of  fundamental  random  variable  w  ,  and  does 

n 

not  involve  a  or  b,  all  scale  factors  involving  constants  a  and  b  in  (7)  have 
disappeared  in  output  z  in  (32),  under  hypothesis  H  .  This  means  that  the 
false  alarm  probability  P  in  (6)  cannot  depend  on  a,b;  put  another  way, 
the  false  alarm  probability  for  the  log -norma  1 i zer  of  figure  1,  subjected  to 


v*«v wn 


the  class  of  inputs  given  by  (7),  is  the  same  for  all  members  of  the  class, 

regardless  of  the  values  of  a  and  b.  Since  the  sample  mean  and  sample 

standard  deviation  In  (32)  still  depend  on  N,  as  seen  by  reference  to  (24), 

the  false  alarm  probability  will  necessarily  be  a  function  of  N,  as  well  as 

depend  on  the  particular  probability  density  function  of  independent 

identically-distributed  normalized  random  variables  v  and  fv  However, 

o  v  n* 

In  general,  there  will  be  no  need  to  investigate  the  false  alarm  probability 
for  the  general  Weibull  class  in  (10),  but  instead  we  can  confine  attention 
to  the  fundamental  exponential  probability  density  function  of  wn  as  given 
by  (8).  Of  course,  vn  must  then  be  the  normalized  random  variable,  as 
given  by  (16),  (17),  (21),  and  (22).  More  details  on  the  statistics  of 
Weibull  variates  and  their  logarithmically-distorted  counterparts  are  given 
in  appendix  A. 

A  similar  statement  can  be  made  with  regard  to  the  fundamental 
log-normal  exceedance  distribution  function  given  by  (11).  In  fact,  the 
logarithmical ly  transformed  input,  (16),  to  the  normalizer  has  exceedance 
distribution  function 

Q^(u)  =  Prob(7  >  u)  =  Prob(jln  w  >  u)  =  Prob(w  >  exp(u))  = 

=  0,  (exp(u))  =  $(-u)  for  all  u  .  (33) 

But  this  is  the  exceedance  distribution  function  of  a  zero-mean 
unit-variance  Gaussian  random  variable.  Thus,  v  of  (16)  is  already  a 
normalized  random  variable,  and  Gaussian  at  that.  Theref ore,  dec i s ion 


variable  z  in  (32)  involves  a  collection  of  N+l  independent  identically- 
distributed  zero-mean  unit-variance  random  variables  vq  and  ^v^  1 . 

Again,  the  false  alarm  probability  can  only  depend  on  N,  and  not  on  scale 
factors  a  and  b  in  (13)  and  (14),  Of  course,  the  detection  probability  (6), 
as  applied  to  (28),  will  depend  additionally  on  parameters  d  and  r  in  (29) 
and  (30). 

In  summary,  the  log-normal i zer  in  figure  1  will  possess  constant  false 
alarm  rate  properties,  that  is,  the  same  false  alarm  probability  for  all  the 
members  of  the  class  of  random  variables  generated  according  to  (7), 
regardless  of  the  values  of  a  and  b. 


PERFORMANCE  FOR  LOG-NORMAL  INPUT  X 


The  problem  of  Interest  in  this  section  Is  the  evaluation  of  detection 
probability  (6)  for  the  decision  variable  z  given  by  (28),  when  normalized 
random  variables  {vj  1  are  independent  Identically-distributed  zero-mean 
unit-variance  Gaussian  random  variables;  this  is  the  case  discussed  in  (33) 
et  seq.  Although  the  probability  density  function  of  normalized  random 
variable  vq  is  arbitrary,  we  will  also  take  it  here  to  be  zero-mean 
unit-variance  Gaussian.  Reference  to  (15),  (17),  and  (25)  reveals  that  this 
is  tantamount  to  assuming  that  the  normalizer  input  {y^  in  figure  1  is 
Gaussian  with  arbitrary  mean  and  variance,  while  yQ  is  also  Gaussian  with 
different  arbitrary  mean  and  variance.  All  these  arbitrary  parameters  are 
collected  together  in  (28)  in  the  parameters  d  and  r,  according  to  (30). 

This  situation  is  also  equivalent  to  assuming  log-normal  excitations  at  the 
input  of  the  log-normal i zer  of  figure  1. 

From  (6)  and  (28),  since  a  >  0, 

PQ  =  Prob(z  >  T)  =  Probed  ♦  rvQ  -  C(v)  >  T  3(v))  ,  (34) 

where  p(v)  and  o(v)  are  given  by  (24).  It  is  shown  in  appendix  B  that  G(v) 
and  3(v)  are  stat ist ica 1 ly  independent,  with  probability  density  functions 
given  by  (B-16)  and  (B-20),  respectively,  (setting  p  =  0,  o  =  1)  as 
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Pp(u) 


(2»/N)  1/2  exp(-  j  u2)  for  all  u 


P  Cu) 

O 


2  mm  u^exof-mu2! 


for  u  >  0;  m  -  *y-  . 


Now  the  quantity  d  +  rvQ  -  y(v)  in  (34)  is  a  Gaussian  random  variable 

2  2 

with  mean  d  and  variance  r  +  1/N  *  aH;  see  (B-16)  or  (35).  Considering  a 
fixed  for  the  moment,  the  conditional  detection  probability  in  (34)  is  then 

r$  L  2on  J  V  N  / 


upon  use  of  (12).  Averaging  this  result  over  the  probability  density 
function  (36)  of  $,  we  have  the  unconditional  detection  probability 


06 

=  J  du  ^ 


<1  -  Tu\  2  mm  u^expf-mu2! 

/  n««) 


where  N  >  2  and 


w"  c  exp( -w  /21 


r{¥] 


J(d;  -  t-  w)  , 


i  _  d  d 


d'  5  - 
r  a 


N  V r  +•  l/N 


°N  ^N-l’l/r2  +  1/N 


The  fundamental  parameters  upon  which  PQ  depends  are 

N,  number  of  normalizer  samples; 

T,  threshold  at  normalizer  output; 
d,  deflection  criterion  (30); 

r,  scaling  (30)  .  (40) 

However,  they  show  up,  in  integral  result  (28),  collapsed  into  the  three 

variables  N,  d' ,  T 1 . 

r  r 

For  signal  absent,  we  have  d=0  and  r=l ,  as  noted  in  (31).  Then  (38)  and 
(39)  reduce  to  the  false  alarm  probability 


which  depends  only  on  N  and  T.  Thus,  given  a  particular  number  N  of 
normalizer  samples,  threshold  T  can  be  selected  to  realize  a  specified  value 
of  false  alarm  probability  P^.  This  applies  for  the  complete  class  of 
log-normal  inputs,  (13)  or  (14),  into  the  log-norma 1 i zer  in  figure  1,  and 
can  be  achieved  without  knowledge  of  a  or  b. 


As  N  ->  <*>  ,  we  have 


°  *  1 .  °N  *  r 


as  N  -*a6  , 


(42) 


1  1 
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giving  from  (37), 

P0  =  §(17J:).  PF  -  f(-T)  for  N  =  «6.  (43) 

This  yields 

f  (PQ)  -  f  [d  +  |( PF)J  for  N  =  (44) 

M 

where  $  is  the  inverse  ^-function;  this  last  result  is  useful  for  plotting 
receiver  operating  characteristics  on  norma  1 -probabi 1 ity  paper.  It 
illustrates  that  for  N=  oo,  those  curves  are  straight  lines  with  slope  l/r 
and  offset  d/r  at  PF=.5. 

The  actual  numerical  evaluations  of  false  alarm  probability  (41)  and 
detection  probability  (38)  are  undertaken  in  appendix  C.  The  inputs  to  the 
functions  considered  there  are  the  3  parameters  N,T^,d^  as  given  by 
(39).  rather  than  the  4  fundamental  parameters  N,T,d,r  listed  in  (40).  This 
is  no  limitation,  since  for  any  given  values  of  N.T.d ,r,  the  quantities 
N.T^.d^,  can  be  easily  calculated  via  (39)  and  used  as  inputs  to  the 
procedures  in  appendix  C. 


I’Ml 


m 
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PERFORMANCE  FOR  WEIBULL  INPUT  * 


Here  we  want  to  evaluate  detection  probability  (6)  for  the  decision 


variable  z  given  by  (28),  when  normalized  random  variables  fv  ]  are 


independent  identically-distributed  zero-mean  unit-variance  log-distorted 


Weibull  variates.  Before  we  do  that,  we  observe  that  detection  probability 


(6)  can  be  expressed  generally  as 


/d  +  rv  -  u(v)  \ 

Pn  =  Prob(  z  >  T)  =  Q  ( T)  =  Prob  - 2 -  >  T  = 

0  z  V  S(*>  / 


,  prob^  >  .  Jdu  \fLr4)  ptlu)  • 


where  we  used  (28)  and  defined  random  variable 


t  =  tl( v)  +  T  S(v) 


in  terms  of  the  sample  quantities  in  (24).  The  separation  of  functions  in 


the  last  form  in  (45)  is  due  to  the  fact  that  random  variables  vq  and  t 


are  statistically  independent  of  each  other.  When  signal  is  absent,  then 


d=0,  -r=l  according  to  (31),  and  (45)  reduces  to  false  alarm  probability 


du  Qu  (u)  p  (u) 
o 


In  the  special  case  where  N  -»«®,  that  is,  a  very  large  number  of  samples 


used  in  the  normalizer,  the  sample  quantities  in  (24)  approach  the  true 
mean  0  and  standard  deviation  1  of  normalized  random  variables  {vp3 ,  and  t 
tends  to  the  constant  T.  Then  (45)  and  (47)  reduce  to 

PD  =  Qv  rr^)  ’  PF  =  Qv  (T)  for  N  =  (48) 

o  '  •  o 

This  limiting  case  can  be  used  as  a  comparison  with  practical  cases  where  N 
is  large,  but  not  infinite. 


We  now  specialize  the  above  general  results  to  the  case  of  log-distorted 
Weibull  variates  {vni 1 .  Although  the  probability  density  function  of 
normalized  random  variable  vq  is  arbitrary,  we  will  also  take  it  here  to 
be  a  normalized  log-distorted  Weibull  variate.  In  this  case,  the  exceedance 
distribution  function  of  random  variable  vq  is  given  by  (A-17)  as 


Qv(u)  =  exp 


for  al 1  u  , 


(49) 


where  y=. 57721  is  Euler's  constant.  Thus,  the  detection  and  false  alarm 
probabilities  for  N  =  °o ,  as  given  generally  by  (48),  are  immediately 
available  upon  use  of  (49). 


The  probability  density  function  of  random  variable  t  defined  in  (46)  is 
a  much  more  difficult  task  for  finite  N.  To  make  any  analytic  progress,  we 
have  had  to  assume  that  t  is  Gaussian;  this  can  be  expected  to  be  a  fair 
approximat ion  if  the  number  of  normalizer  samples  N  entering  the  sample 


20 


quantities  in  (24)  is  large,  according  to  the  central  limit  theorem. 
However,  we  can  anticipate  that  the  approach  of  random  variable  t  to 
normality  will  be  faster  near  its  mean,  but  considerably  slower  on  the 
tails.  This  can  lead  to  a  significant  bias  in  the  calculation  of  small 
false  alarm  probabilities. 


Thus,  our  assumption  is  that  the  probability  density  function  of  t  in 
(46)  is  given  by 

Pt(u)  =  [2»  oW'2  expf-  *  ■ 

1  L  2  a‘(t) 


(50) 


Combining  (49)  and  (50)  in  (45),  the  detection  probability  is  given 
(approximately)  by 


2  1 

ml 

(t)  - 


«  (2ir) 


|i.  .*(«]■  .,,[-1*^ 

L  J  L  2  <T(t 

1/2  J  dx  exp^"-  -exp^  +•  h?  x^ 


(51) 


where  constants 


IT  u(  t)  -  d 

r 


=  — it  mi 
2  r 


(52) 


The  false  alarm  probability  follows  upon  setting  d  =  0,  r  =  1  in  (52). 


The  fundamental  parameters  in  integral  result  (51)  are  d  and  r,  along 


with  mean  y(t)  and  standard  deviation  o(t)  of  random  variable  t.  The 
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complexity  of  random  variable  t,  defined  by  (46)  and  (24),  precludes  us  from 
evaluating  mean  y(t)  and  standard  deviation  a(t)  exactly.  However,  numerous 
simulations,  each  consisting  of  100,000  trials,  enabled  us  to  extract  the 
following  rather  accurate  rules  of  thumb  for  the  statistics  of  t. 


9 


First  of  all,  for  general  definition  (46),  we  have  mean 

P(t)  =  y{C(v)J  +  T  p{S(v)}  ,  (53) 

and  variance 

a2(t)  =  a2{0(vj}  +  T2  «2{S(V)}  2T  of{J(  V  jj  o{S(  V )]  p  ,  (54) 

where  p  is  the  normalized  correlation  coefficient  between  tl(v)  and  3(v). 

The  simulation  results  alluded  to  above,  for  normalized  random 
variables  fvn}  in  (24)  being  log-distorted  Weibull  variates,  are  given  by 

v{v(vj}  =  0  ufc(v)]  M  -  y 

o2fu(v)]  .  1  ^(V)l  . 

P  2  -.55  .  (55) 

Observe  the  large  value  of  p,  in  contrast  with  the  earlier  case  of  a 
Gaussian  input  to  the  normal i zer,  where  the  sample  mean  and  sample  standard 
deviation  were  not  only  uncorrelated  but  in  fact  independent;  see  appendix  8. 


r. 

< 
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The  quantities  in  (55)  depend  solely  on  N;  when  used  in  (53)  and  (54), 
it  is  seen  that  y(t)  and  o(t)  in  (54)  depend  on  both  N  and  threshold  T. 
Thus,  the  totality  of  fundamental  parameters  of  relevance  in  detection 
probability  (51)  is  d.r.N.T,  just  as  for  the  Gaussian  case  in  (40). 

Analytic  evaluation  of  Integral  (51)  is  impossible;  accordingly,  numerical 
Integration  was  employed. 

An  alternative  approximation  to  the  statistics  of  random  variable  t  of 
(46)  is  undertaken  in  appendix  0.  Namely,  for  large  threshold  values  T, 
where  random  variable  t  is  dominated  by  sample  standard  deviation  9(v),  it 
might  be  thought  that  a  X-approximation  would  have  wider  applicability  than 
a  Guassian  one.  This  is  indeed  true,  as  will  be  demonstrated  by  the 
numerical  results  to  follow;  however,  for  small  false  alarm  probabilities, 
the  x-approximation  also  falls  short. 


■>/lVrvvA'AV.V.v‘v^ 


?  3 
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GRAPHICAL  RESULTS 

In  this  section,  we  will  present  graphical  results  for  the  exceedance 

distribution  functions  and  receiver  operating  characteristics  for  both  the 

Gaussian  and  the  log-distorted  Weibull  inputs  to  the  normalizer  in  figure  1. 

This  corresponds,  respectively,  to  log-normal  and  Weibull  inputs  (that  is, 

detector  outputs)  to  the  logarithmic  device  in  figure  1.  The  theoretical 

results  of  the  previous  two  sections  are  augmented  by  simulation  results, 

23 

each  based  upon  2  =  8.4  million  trials  of  decision  variable  (normalizer 

output)  z  given  by  (28)  and  (24). 

The  scaling  parameter  r,  that  is,  the  ratio  of  standard  deviations  in 
(30),  is  taken  at  1  for  all  these  results,  in  order  to  keep  the  number  of 
plots  at  a  reasonable  level.  The  deflection  parameter  d  in  (30)  is  varied 
from  0  to  values  large  enough  to  sweep  out  the  important  range  of  detection 
and  false  alarm  probabi 1 ities  of  interest.  The  number  of  normalizer 
samples,  N,  is  taken  at  the  values  « ,  64,  32,  16,  which  appears  to  cover 
the  most  important  range  of  practical  use.  Threshold  value  T  in  exceedance 
probabilities  (6)  is  allowed  to  vary  widely,  so  that  the  full  range  of 
detection  and  false  alarm  probabilities  can  be  observed. 
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GAUSSIAN  INPUT  TO  NORMALIZER 


The  exceedance  distribution  function  (EOF),  defined  in  (6),  for  N=«o  is 
plotted  in  figure  2  versus  threshold  T,  for  deflection  parameter  d  taking  on 
values 

d  =■  0(1)9  *  0.1, 2, 3. 4, 5, 6, 7, 8,9  .  (56) 

The  arrow  on  the  figure  indicates  the  direction  of  increasing  d;  thus  d=0, 
which  is  the  false  alarm  probability,  corresponds  to  the  curve  at  the  lower 
left.  The  results  in  figure  2  are  based  on  (43);  since  the  ordinate  is 
according  to  a  normal  probability  rule,  these  curves  are  perfectly  straight 
lines.  Exceedance  probabilities  ranging  from  IE-6  to  .999999  are  covered 
when  threshold  T  is  varied  over  the  range  -5,5. 

When  N  takes  on  the  values  64,  32,  16,  the  corresponding  results  are 
displayed  in  figures  3,  4,  5,  respectively.  These  graphs  were  obtained  from 
( 38 ) -( 4 1 ) .  implemented  by  the  procedures  in  appendix  C.  The  exceedance 
distribution  function  for  N=64  in  figure  3  is  fairly  close  to  that  for  N=» 

In  figure  2;  however,  by  the  time  the  number  of  normal izer  samples  N  has 
decreased  to  16  in  figure  5,  significant  curvature  has  developed  in  the 


Exceedance  Probability 
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.399999 
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Superposed  in  figure  5  are  11  simulation  results  for  d=0( 1)10,  each 


based  upon  8.4  million  Independent  trials.  Due  to  the  large  number  of 


trials,  the  theoretical  and  simulation  results  are  indistinguishable,  except 


near  the  extremes  of  probability  IE-6  and  .999999,  where  the  jagged 


character  of  simulation  results  is  manifested.  This  close  agreement  of 


results  not  only  confirms  the  theoretical  analysis  but  also  lends  credence 


to  the  use  of  simulation  for  the  estimation  of  probabilities  out  on  the 


tails  of  the  distribution,  provided  that  enough  trials  are  conducted. 


Figures  3,  4,  5  furnish  information  which  enables  the  selection  of  the 


required  threshold  T  to  realize  a  specified  false  alarm  probability  for 


N  =  64,  32,  16,  respectively.  For  example,  figure  5  with  d  =  0  indicates 


that  to  realize  a  false  alarm  probability  of  IE-5  for  N  =  16,  threshold  T  in 


(6)  must  be  chosen  as  6.3. 


When  threshold  T  is  eliminated,  and  the  detection  probability  plotted 


versus  the  false  alarm  probability,  we  obtain  the  receiver  operating 


characteristics  (ROC).  The  result  for  N=o«  is  given  in  figure  6,  where  both 


the  abscissa  and  ordinate  are  plotted  according  to  a  normal  probability 


scale.  Deflection  parameter  d  varies  over  the  range 


d  =  0< .5)7.5  =  0,. 5, 1,1. 5 . 7,7.5  . 


The  arrow  again  points  in  the  direction  of  increasing  d;  thus  d=0  is  the 


curve  on  the  lower  right.  These  curves  are  precisely  the  straight  lines 


indicated  by  (44)  with  r=l. 
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Corresponding  receiver  operating  characteristics  for  N  =  64,  32,  16  are 
presented  in  figures  7,  8,  9.  The  detection  and  false  alarm  probabi 1 ities 
both  range  down  to  IE-6,  while  the  upper  limits  have  been  truncated  at  .99 
and  .5,  respectively.  Values  beyond  these  limits  can  be  obtained  from  the 
earlier  figures  2  through  5. 

Superposed  in  figure  9  are  ten  simulation  results  for  d  =  1(1)10. 

Again,  except  for  the  small  probability  regions  like  Pp  <  IE-5,  the 
theoretical  and  simulation  results  are  indistinguishable  and  overlay  each 
other.  It  will  be  noticed  that  a  characteristic  wiggle  in  the  receiver 
operating  characteristics  is  duplicated  for  every  simulation  result,  at  a 
constant  value  of  false  alarm  probability;  for  example,  see  the  triangular 
bump  in  all  10  simulation  results  at  Pp  a  IE-6.  The  reason  for  this 
behavior  is  that  when  random  variable  z  in  (28)  was  simulated,  the  random 
numbers  employed  in  (24)  to  generate  y  and  a  were  not  changed  when  different 
d  values  were  considered  in  (28).  The  reason  for  this  deliberate  choice  was 
economy  of  computer  execution  time;  that  is,  the  time-consuming  task  of 
computation  of  (24)  was  done  once  for  each  trial,  and  used  in  (28)  for  all 
of  the  d  values  of  interest.  This  repeated  use  of  the  same  y,o  values  for 
different  d  values  gives  a  persistent  systematic  perturbation  to  the 
estimated  receiver  operating  characteristics  at  a  fixed  false  alarm 
probability.  However,  for  8.4  million  trials,  this  bias  is  small,  even  for 
the  rare  events  with  probabilities  greater  than  IE-6,  and  was  deemed 
acceptable  in  light  of  the  greatly  increased  computer  time  required  for  the 
alternative  approach  of  regeneration  of  y  and  a. 


Figure  7.  Probability  of  False. Alarm 
ROC  for  N  =  64;  Gaussian 
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Probability  of  Detection 


Probability  of  Detect 
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As  an  example  of  the  use  of  figures  6  through  9,  the  values  of  deflection 
d  required  to  realize  =  IE-5  and  PQ  =  .5  are 

d  =  4.3,  4.7,  5.1,  6.2  for  N  =  *6  ,  64,  32,  16  ,  (58) 

respectively.  The  cost  of  reducing  N  from  to  16  is  that  d  must  be 
increased  by  the  factor  6. 2/4.3  =  1.44;  whether  this  is  tolerable  depends  on 
the  application.  The  relation  of  deflection  parameter  d  to  any  system  input 
slgnal-to-noise  ratio  depends  on  the  particular  processor  form  preceding  the 
logarithmic  device  in  figure  1,  and  must  be  left  to  the  user  and  his 
particular  application. 

LOG-WE IBULl  INPUT  TO  NORMALIZER 

When  the  input  to  the  normalizer  is  a  log-distorted  Weibull  variate,  the 
performance  is  markedly  different.  The  exceedance  distribution  function  for 
N=  <*>  is  displayed  in  figure  10  and  has  a  significant  curvature  when  plotted 
on  normal  probability  paper;  these  results  are  based  upon  the  use  of  (48) 
and  (49).  The  use  of  the  notation  'Extreme1  is  explained  in  (A-7)  et  seq. 

When  N  is  decreased  to  64,  the  corresponding  exceedance  distribution 
functions  are  given  in  figure  11.  Due  to  the  questionable  assumptions 
required  in  the  theoretical  analysis  of  this  case  and  used  in  (50)  et  seq., 
simulation  results  were  also  superposed  for  the  values  d  -  0 ( .5)5. 

Agreement  in  the  mid-range  of  probabilities  is  excellent.  At  the  low  end  of 
the  probability  range,  near  IE-6,  the  simulation  results  indicate  a 
systemat lea  1 ly  lower  exceedance  probability  than  predicted  by  theory. 
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This  erratic  trend  of  the  theoretical  approximation  is  continued  and 
accented  in  figure  12  for  N=32  and  in  figure  13  for  N=16.  In  fact,  in  the 
latter  case,  for  threshold  T=5,  the  simulation  indicates  exceedance 
probabilities  for  d=0  that  are  more  than  2  orders  of  magnitude  smaller  than 
the  theory  predicts;  see  bottom  right  of  figure  13.  The  discrepancies  at 
the  high  end  of  probabi 1 ities  are  also  considerable,  as  seen  at  the  top  left 
of  the  figure. 

Also  added  to  this  particular  figure  is  the  result  of  using  the 
X-approximation  for  random  variable  t  of  (46),  as  detailed  in  appendix  0. 
Although  the  improvement  in  probability  values  is  over  an  order  of  magnitude, 
there  is  still  another  order  of  magnitude  error  left  in  this  alternative 
approximate  approach.  The  reason  for  the  difficulty  in  the  theoretical 
analysis  is  two-fold:  (1)  values  of  N  like  16  or  32  are  not  large  enough 
for  the  central  limit  theorem  to  have  developed  substantial  accuracy  on  the 
tails;  (2)  the  probability  density  function  of  a  log-distorted  Weibull 
variate,  as  given  by  (A-7),  is  distinctly  non-Gaussian  on  the  tails.  The 
decay  of  (A-7)  on  the  positive  tail  is  much  faster  than  Gaussian,  while  that 
on  the  negative  tail  is  slower,  being  only  exponential. 
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The  receiver  operating  characteristics  for  N=«o  are  given  in  figure  14, 
while  those  for  N  =  64, *32,  16  are  given  In  figures  15  through  17, 
respectively.  The  discrepancy  between  theory  and  simulation  becomes 
progressively  larger  as  N  decreases,  reaching  the  point  in  figure  17  where 
the  theory  is  entirely  invalid  for  false  alarm  probabilities  less  than 
approximately  .001.  The  reason  for  the  severe  dip  of  the  theoretical  curves 
to  the  left  of  each  figure  is  the  inadequacy  of  the  false  alarm  probability 
approximation,  it  being  much  too  large  for  the  larger  threshold  values;  see 
bottom  right  of  figure  13.  On  the  other  hand,  the  simulation  results  in 
these  figures  are  all  based  on  8.4  million  independent  trials,  making  them 
trustworthy  well  down  near  the  1 E  —6  level  of  probability  plotted  here. 

As  an  example  of  the  use  of  figures  14  through  17,  the  values  of  d 
required  to  realize  probabi 1 i ties  Pf  «  1 E -5  and  PQ  =  .5  are 

d  =  2.2,  2.6,  2.9,  3.8  for  N  =  «©  ,  64,  32,  16  ,  (  59) 

respectively.  The  latter  three  values  are  extracted  from  the  simulation 
results  in  figures  15  through  17.  Oirect  comparison  of  the  absolute  levels 
in  (59)  with  the  correspond i ng  Gaussian  results  in  (58)  is  not  valid, 
because  the  shapes  of  the  input  probability  density  functions  in  the  two 
cases  are  markedly  different  and  are  more  important  than  the  deflection 
criterion,  defined  by  (30). 
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Figure  14.  Probability  of  False  Plarm 
ROC  for  N  * » ;  Extreme 
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Figure  15.  Probability  of  False  Rlarm 
ROC  for  N  =  64;  Extreme 
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CONCLUSION 

The  performance  of  cne  normalizer  with  a  Gaussian  input  is  capable  of 
exact  analysis  in  terms  of  integrals  which  are  readily  evaluated  via 
recursions.  The  main  reason  that  this  fortuitous  situation  obtains  is  the 
statistical  independence  of  the  sample  mean  and  sample  standard  deviation 
for  Gaussian  random  variables.  However,  for  other  inputs  to  the  normalizer, 
these  sample  statistics  are  highly  correlated  with  each  other  and  create  an 
untractable  analysis  problem. 

An  acceptable  alternative  in  this  latter  case  is  simulation  with  a  large 
number  of  trials.  Here  8.4  million  trials  were  employed,  which  allowed  for 
estimation  of  tail  probabi 1 ities  in  the  1 E -6  range.  If  the  false  alarm 
probability  could  be  evaluated  theoretically,  then  simulation  would  only 
need  to  be  conducted  for  the  detection  probability  P^.  And  if  PQ  were 
of  interest  only  in  the  range  (.5, .99)  say,  then  as  few  as  10,000  trials 
would  suffice  for  a  decent  estimate.  However,  it  appears  that,  in  general, 
even  the  analysis  for  the  false  alarm  probability  involves  some  unmanageable 


statistical  relations. 
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APPENDIX  A.  WEIBULL  VARIATES 


The  exceedance  distribution  function  of  a  general  Weibull  random 
variable  x  is  given  by  C?*P-^21 


Qx(u)  =  Prob(x  >  u)  =  exp 


t  ftii 


for  u>0;a>0,  b>0.  (A-l 


The  corresponding  probability  density  function  of  x  is 


,  >L-’ 
PX(U)  ■  -Qx(u)  = 

ab  b 


exp 


1/b 


for  u  >  0  . 


( A-2  ] 


The  w-th  moment  of  x  is 


t”  =  J'  du  uu  px( u )  =  av  f(l  +-  bv)  for  v  >  -1/b  . 


( A-3) 


The  characteristic  function  of  x  is  not  available  in  closed  form,  for 


general  b;  however 


exp(ifx)  =  (1  -  iTa)-1  for  b  =  1  .  ( A -4 ) 

The  normalized  cumulants  of  x,  for  general  b,  are  independent  of  a;  however, 
they  do  not  approach  zero  as  either  b  ->  0  or  b  -» °o  .  Therefore  x  does  not 
tend  to  Gaussian  as  the  shape  parameter  b  is  changed. 


TR  9015 


LOG-OISTORTEO  WEI8ULL  VARIATE 

As  indicated  in  (1)  and  (15),  we  are  interested  in  the  log-distorted 
random  variable 

y  =  In  x  ,  ( A— 5 ) 

where  x  is  a  Weibull  variate  with  probability  density  function  (A-2).  The 
exceedance  distribution  function  of  y  is 

Qy(u)  =  Prob(y  >  u)  *  Prob(ln  x  >  u)  =  Prob(x  >  exp(u))  = 

for  all  u  ,  (A -6) 


L  a  . 


where  (A-l)  was  employed.  The  corresponding  probability  density  function  of 
random  variable  y  is 


Py(u)  =  -Qy(u)  = 


b  a 


1/b 


exp 


u 

_b 


exp( u/b) 
1/b 


for  all  u 


(A-7) 


which  is  a  form  of  the  probability  density  function  for  extreme  values;  see 
[4;  (14.65)].  We  will  refer  to  (A-7)  as  an  extreme  value  probability 
density  function  here. 
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The  characteristic  function  of  tandom  variable  y  is 

fy(i?)  =  exp(ify)  =  exp(if!n  x)  =  x1*  = 

=  air  T(1  +  Ifb)  ,  ( A —8 ) 

the  last  step  by  use  of  (A-3);  this  is  a  generalization  of  [4;  page  344, 
exercise  14.4].  The  actual  numerical  evaluation  of  the  characteristic 
function,  in  £ A— 8 )  for  real  f  is  best  accomplished  by  employing  (A-7): 


fy( i?)  =  exp(ijy)  =  j  du  exp(i?u)  py(u)  = 

■Ats  f  «"«*<«»>  «» [» -  “Hf1]  •  <*-»> 

b  a  J  L  a  J 

This  can  be  efficiently  and  accurately  evaluated  by  use  of  a  fast  Fourier 
transform;  the  integrand  decays  very  rapidly  as  u  -»  ±  . 

In  anticipation  of  getting  the  cumulants  of  random  variable  y,  we  have 
from  ( A-8) , 

In  fy(i?)  =  i?  In  a  +  In  P(1  +  iTb)  .  (A-10) 

Now  from  [5:  (6.1.33)  and  section  23.2]  , 
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In  P(1  +  z)  *  -yz  ♦  (-l)n  }{n)  zn/n  , 

n-2 

where  y  =  .57721  is  Euler's  constant  and 

OO 


(A-11) 


( A— 1 2) 


k=l 


In  particular,  ^(2)  =  *^6. 


There  then  follows,  from  (A-10)  and  (A-11),  the  cumulants  of  random 
variable  y  as 


In  a  -  by 


for  n=l 


*y<">  ' 


( A— 1 3) 


(-1)"  ?(n)(n-l)!  bn  for  n  >  2 

2  2/ 

In  particular,  the  variance  of  y  is  Xy(2)  =  b  w/6.  For  n  >  2,  the 
normalized  cumulant  of  y  is 


%y(") 


Wy(2>] 


n/2 


•  (-  ®)n  . 


( A  —  1  4 ) 


which  is  independent  of  both  a  and  b;  thus  random  variable  y  does  not 
approach  Gaussian  as  a  and/or  b  approach  any  limits  whatsoever.  These 
results  generalize  [4;  page  344,  exercise  14.4], 


A  4 


aVs»‘  . 


.'.V.V 


4 

j 

I 
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NORMALIZED  LOG-DISTORTED  WEIBULL  VARIATE 


If  a=b=l  in  (7),  then  xn  =  wn,  and  it  then  follows  from  (16)  and  (1) 

that  =  In  w^  =  In  xn  =  yn>  In  this  case,  we  can  use  (A-6)  with  a=b=l  to 

obtain  the  exceedance  distribution  function  of  v  as 

n 


Q^(u)  =  exp  [”  -exp(u)]  for  all  u  . 


(A-15) 


Additionally,  there  follows  from  (A-13) 


w(v)  =fyl)  =  -y  =  -.57721  , 
«(v)  =(Xy(2 ))=  ir/yT  . 


(A-16) 


We  are  now  in  position  to  determine  the  exceedance  distribution  function 

of  the  normalized  log-distorted  Weibull  random  variable  v  defined  in 

n 

(17),  namely, 


Qv(u)  =  Prob(v  >  u)  =  Prob(v  >  y(v)  +  a(v)u)  = 

=  Qy(y(v)  +■  o(v)u)  =  exp  r~exp(-Y  +  ~  u)  for  all  u  .  ( A -1 7 ) 

L  Y6 


A  -  b 
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LOG-NORMAL  VARIATES 

For  completeness,  we  list  here  the  v-th  moment  of  log-normal  variate  x 
with  probability  density  function  as  given  by  (14)  and  (12): 

Oft 

X*  =  f  du  u»  Px(u)  - 1  J  as  „• .  = 

-«o  ( A-l  8) 


A-b 
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■  ■  .U  .1.  -t.  JU» 
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APPENDIX  B.  INDEPENDENCE  OF  SAMPLE  MEAN  AND  SAMPLE  VARIANCE 


FOR  GAUSSIAN  RANDOM  VARIABLES 


Let  fxj”  be  independent  identically-distributed  Gaussian  random 


variables  with  mean  and  variance 


—  2  2 

xn  =  v  .  (xn  -  u)  a  a  for  all  n  . 


( B-l  ) 


Define  sample  mean 


( B-2 ) 


and  sample  variance 


(B-3) 


where  scale  factor  g  =  1/N  or  1/(N-1)  typically.  (N  >  2  required.) 


We  have,  in  vector  notation 


m  =  ~  1T  x  =  J[1  1 


.  1  ][x1  x2  . . .  x^]  , 


(B-4) 


v  =  g 


r N  1 

2  xn  *  m2  N  =  g[*T*  ~  q  XT1  lTx]  =  g  X1 
.  n=l 


(B-5) 


where 


Q  -  I  -  1  1  1T 


The  joint  characteristic  function  of  m  and  v  is 


f(f,d)  =  exp(i^m  +  iev)  = 


Now  the  joint  probability  density  function  of  vector  X  is 


N  r> 

p(x) =  Tf  exp  ■ 

n=1  L 


(xn  -  u)‘ 


2  a 


=  (2ir  a2)  exp  p  (XTX  -  2v  1 TX  +  p2nJJ  . 


Therefore 


-N/2  /*  r 

f  (F,e)  =  (2ir  a)  dX  exp  -  — (x'x-2w  IX  +  p  n)  + 

'J  k  2a 


l  2a 


.  9  T  T 

+  i  jj  I  X  +  i«g  X  QX 


Now  we  use  [6;  (8-1 ) ] 


F  (  2ir ) N 

1/2 

1  ,  T  U-1  , 

[ det  M j 

exp 

2  L  M  L 

J 

with  identifications 


1 


M  =  —  I  -  i2eg  0  .  L  = 

a 


-H  i.  i  f 


2  ♦  1  N  1 


8-2 


(B-b) 


=  expfifj  1TX  +  ieg  XTQx]  .  (B-7) 


(  B — 8 ) 


(B-9) 


(8-10) 


(8-1  1 ) 
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Using  the  definition  of  0  in  (B-6),  there  follows 


M  *  -  12«g^  I  +  ^  1  1T  . 


(B-12) 


Now  from  [6;  (21)  and  (22)], 


det  M  =  -—(1  -  i2«J2eg)N"1 

<3 


.  2 

l  -  i2o  eg 


[,  .  ilsaa!  ,  ,t] 


( B-l 3) 


There  follows 


(B-14) 


f ( T. ® )  =  exp  i 


(B-l 5) 


Since  this  joint  characteristic  function  factors,  it  follows  that  sample 
statistics  m  and  v  are  stat i st lea  1 ly  independent.  Also  the  probability 
density  functions  are  obviously 


Pm(>J)  =  ~~  ~  ~  e*P  -  ^L 

V27  o/fiT  L  2 


for  all  u 


( B-l  b) 


Pv(u) 


u  expl?a?c 


r(¥'  (*4 


tor  u  >  0 


(B-l  /) 
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2 

Thus  sample  mean  m  is  Gaussian  with  mean  p  and  variance  o  /N;  while  sample 

2 

variance  v  is  chi-squared  of  N-l  degrees  of  freedom  with  mean  a  (N-l)g  and 
4  2 

variance  2a  (N-l)g  .  For  the  typical  choice  of  gain  g  =  1 / ( N-l ) ,  this 

2  4 

implies  that  v  has  mean  a  and  variance  2a  /(N-l),  and  therefore 


lim  pv(u)  =  4(u  -  a  ) 

N-»«* 


( B — 1 8) 


The  sample  standard  deviation 

s  =  'fT 

has  probability  density  function 


(II)  •  2  u  Pv(u2)  •  ?  a"  for  u  >  0 


where 


2a 


( B  —  1  9 ) 


( B-20) 


0  =  for  ^  =  1 /(N-l ) 


(8-21) 


The  k-th  moment  of  s  is 


oo 


k  f .  •  k  ,  .  f  . 

s  =  1  du  u  p$(u)  =  du 


2  uN~2^»rC-u2/B) 
s'-'  1  N-l 

°  rMs2 


r/N±k-j\  k/; 

ilZOJl  ... 

r(¥) 


.  ( B  -22 ) 


In  particular, 


2  2 
s  =  o 


( B  -23 ) 


8  -4 


TR  8075 


and  [5;  (6.1.47)] 


s  =  a 


'»  iff) 

rm 


/  3\  1/2  r 


i  ♦  -i^V  +  o 

/  n\2 


H) 


’  1/4  1/16 

"I  *"»># 

..r.  -  ♦  o(n-3,] 

N  -  -  J 

11  Q 


>-8 


(<  -  iy 


+  0(N  J) 


as  N  •*  «o  . 


J(N  -  ^)’3_ 
f  - 1)"3~ 


( B-24 : 


85/Bb 
Reverse  Blank 


V.WiViS 


.w 


V*a  j1*  ■  «.-  (.«  ».»  <  ■ » 


vENWSSj 


.niii’.-kLiLst. 
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APPENDIX  C.  PROBABILITY  RECURSIONS  FOR  GAUSSIAN  CASE 


The  detection  and  false  alarm  probabilities  are  given  in  integral  form 


in  (38)  and  (41).  These  integrals  have  already  been  encountered  in 


[1;  appendix  E],  and  evaluated  in  a  recursive  fashion.  We  will  modify 


those  results  somewhat,  in  order  to  better  suit  the  current  forms. 


First,  we  have,  from  (41)  and  (39), 


r  N-2  ,  2... 

o  ,2 


PpCN.T, ) 


(C-l) 


Define 


t1  =  ^1  +  T^2^  ,  where  T^  =  T< 


N2  -  1 


Then  from  [1;  ( E - 1 7 ) ] ,  using  identifications  (that  is,  replacements  from 


there  to  here) 


<  * 


,  K  -*  N-2  , 


there  follows  the  simple  result 


B 


TR  8075 


JL* 

2  £ 


PF(N,Tj)  =  \  -  ;atn(Tj) 


I T!  xi 


bk », 


for  N=2,4,6, 


(C-4) 


k=0 


where 


bo  "  1  * 


bk  =  bk-i  — T 
k  +  2 


for  k  >  1  . 


(C-5) 


A  program  for  this  false  alarm  probability  is  given  in  appendix  E  under  the 
name  FNPf246,  where  t|  is  represented  by  variable  Tp. 


Also, 


N-3 

2 


PF(N’T1)  =  2  “  2  T'lK2ak  xl  f0r  N=3-5’7 - 


(C-6) 


where 


ao  =  1  • 


k=0 


k-i 

a,  =  a.  .  — r —  for  k  >  1  . 
k  k-1  k 


(C-7) 


These  results  are  very  tractable  and  efficient  forms  for  recursive  computer 
evaluation.  A  program  for  (C-6)  is  given  in  appendix  E  under  the  name 
FNPf357,  where  T|  is  represented  by  variable  Tp. 


I 

dT  *  dr  ’  r  *  7  ,V/2  * 

C  +  T'r2) 


K  ->  N-2  . 
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( C— 8) 


("he  curves  in  [1]  are  not  directly  applicable  here  because  they  employed 
the  fundamental  parameter  d^  -»  d*.  which  is  d/(r2  +  1/N)here;  however, 
the  recursions  derived  there  are  immediately  useable.)  We  can  then  use 
[1;  (E-8)]  to  develop  an  expression  for  PQ,  in  terms  of  the  auxiliary 
sequence  (g(k')}  defined  in  [1;  ( E-7 )  ] .  In  particular,  [1;  (E-9)]  yields,  with 


'r  -  (  ♦  T'r2)"’  ■ 


(C-9) 


the  result 


9(01  -  exp  (-  i  d;2xr)  §(o;t;{x7)  ; 


(C-10) 


[1;  ( E-l 3) ]  yields 


—  ■)  j  2 

9(U  -  exp (-d;2/2)  *  (;)  •>'  9(0) J  ; 


(c-n ) 


and  [1 ;  (E-l 2) ]  yields 


g(K)  =  x  h(K)g(K  -  1)  + 


¥  0(K-2>] 


for  K  >  2  , 


( C  —  1  2 ) 


with  definition 


'i/7  rr  r(fn) 


& 


(C-13) 
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Then  we  can  also  use 


h<°)  -  r;  d;(f)  • 


h( K)  =  h(K-2)  —■  for  K  >  2  . 


Finally,  Pq  =  Po(N,Tr,dr)  is  given  by  [1;  (E-8)]  as 


( C— 1 4 ) 


1(0  -  g(K) 

K=0 
K  even 


for  N=3 , 5,7, 


( C-l 5) 


Pn?  “  g(K)  for  N-4,6,8, 
K=1 
K  odd 


where  P02  is  the  value  of  detection  probability  PQ  for  N  =  2.  Observe 
that  the  input  parameters  to  PQ  are  N.T^.d^,  rather  than  the  four 
fundamental  parameters  in  (40);  that  is,  N,T,d,r  are  collapsed  into 
N.T^.d^  according  to  (39).  Programs  for  ( C - 1 5 )  are  furnished  in 
appendix  E  under  the  names  FNPd357  and  FNPd246,  respectively,  where  and 


d^  are  represented  by  variables  Tp  and  Dp. 


The  quantity  PQ2  in  (C-15)  is  evaluated  according  to  the  method  in 
[1;  appendix  F);  an  error  tolerance  and  maximum  number  of  terms  must  also  be 
specified  to  terminate  the  infinite  sum  given  by  [1;  (F-2)]. 


«  *  ,  <r  *r  * 


.  •  ,  ■  /•  «N  A  -  .'i 


TR  8075 


APPENDIX  D.  ^-APPROXIMATION  FOR  RANDOM  VARIABLE  t 


Suppose  we  assume  that  the  random  variable  t  in  (46)  is  a  multiple  of  a 


X-variate  with  K  degrees  of  freedom;  then  its  probability  density  function 


is  [7;  pages  5-7  for  v  =  1/2] 


Pt(u)  = 


uK  1  expj-u2  /  (2A2) 


AK  22 


1  r(l) 


for  u  >  0  . 


Then  the  w-th  moment  of  random  variable  t  is 


tw  =  j'  du  u”  pt( 


2v/2A  u  rm 


and  in  particular 


=  A 


Wt)  7, 


2  2 
V  =  A^K 


Then  the  ratio 


fr§ 


1  -3 

77-  *  0(k  J) 


.v.v.v, 


(0-1) 


(0-2) 


( D-3) 


(0-4) 
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where  the  last  result  uses  the  development  In  (B-24)  with  N  replaced  by  K+l . 


Given  a  value  for  ratio  R  on  the  left  side  of  (D-4),  K  can  be  solved  for 
uniquely,  since  the  ratio  involving  gamma  functions  increases  monotonical ly 
from  0  to  1  as  K  goes  from  0  to  +«a.  In  fact,  to  a  good  approximation  for 
large  K,  the  last  part  of  (0-4)  gives 


K  3  4( 1 -R )  "  8  • 


(0-5) 


Here  we  are  allowing  K  in  probability  density  function  in  (D-l)  to  be 
arbitrary,  that  is,  not  limited  to  integer  values.  Then  we  can  solve  for 
the  required  value  of  A  according  to  (0-3),  as  A2  =  t2A.  This 
procedure  fits  the  assumed  probability  density  function  form  in  (D-l)  to 
specified  values  of  the  first  two  moments  of  t  given  by  (0-3),  as  given  by 
simulation  results  (53)— (55) . 


If  we  now  employ  the  "^-approximate  probability  density  function  for  t 
given  by  (D-l)  in  detection  probability  result  (45),  along  with  (49),  we 
obtain 


P 


0 


l 


du  exp 


-exp  Mr  ♦  -p. 


K-l 

u  exp 


'  K 
=  22 

where  constants 


exp 


[- 


x2/2 


exp(  h 


(0-6) 


h 


1 


ir  d 


(0-7) 


the  numerical  evaluation  of  (0-6)  was  undertaken  for  N  -  16,  and  is 
discussed  in  the  Graphical  Results  section  of  this  report  . 


0  -2 
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APPENOIX  E.  PROGRAMS 

In  this  appendix,  four  programs  are  listed.  They  are  written  in  BASIC 
for  the  Hewlett-Packard  9000  Model  520  Desk  Top  Computer.  Their  titles  are 

EOF  -  Gaussian, 

EOF  -  Extreme, 

Simulation-Extreme, 

Plot-Simulation. 

The  first  one  computes  the  exceedance  distribution  function  for  a  Gaussian 
input  to  the  normalizer  of  figure  1,  for  N=16  (line  10)  and  for  d  =  0(1)12 
(lines  960-970).  This  program  is  heavily  based  on  the  results  of  appendix  C. 

The  second  program  computes  the  approximate  exceedance  distribution 
function  for  a  log-distorted  Weibull  input  to  the  normalizer  of  figure  1, 
for  N= 1 6  (line  10),  r=l  (line  20),  and  d  =  0 ( .75)7.5  (line  1070).  It  is 
based  on  numerical  integration  of  (51)  via  Simpson's  rule. 

The  third  program  simulates  the  normalizer  output  (28)  and  (24)  for  a 

log -distorted  Weibull  input,  for  N-16  (line  10),  d  =  0 ( .75)7.5  (line  20), 

23 

r=  1  (line  30),  and  2  =  8.4  million  trials  (line  40).  The  range  of 
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values  In  z  Is  (-15,5),  which  is  divided  into  1000  bins;  see  lines  60,  90, 
100.  The  resultant  histogram  is  then  summed  on  the  upper  tail  to  yield  the 
exceedance  distribution  function.  The  fourth  program  plots  these  simulation 
results  for  the  exceedance  distribution  function  vs  threshold  T. 
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Table  E-l.  EDF  -  Gaussian 


10  Ns* 1 6  !  NUMBER  OF  SAMPLES 

20  X 1 *-7  !  THRESHOLD 

30  X2*7  !  LIMITS 

40  DIM  A*C303,B*C30] 

30  DIM  X1abel*<l:30>,Ylabe1*<ls30) 

60  DIM  Xcoord<lj30),Ycoord<l830> 

70  DIM  Xgr id<l»30>,Ygrid<l:30> 

80  DOUBLE  Lx,Ly,Nx,Ny, I, Ns, It 

90  ! 

100  A**" Threshold" 

110  B**"E:<ceedance  Probability" 

120  • 

130  Lx*  1 3 

140  REDIM  XI abel S< 1 : Lx) , Xcoordv 1 : Lx > 

130  DATA  -7, -6, -3, -4, -3, -2, -1,0, 1,2, 3, 4, 5, 6, 7 
160  READ  X 1 abe 1  *  <  *  > 

170  DATA  -7, -6, -3, -4, -3, -2, -1,0, 1,2, 3, 4, 3, 6, 7 
130  READ  Xcoord<  *> 

190  ' 

200  Ly-27 

210  REDIM  Y1 abe I *< l : Ly > , Ycoord' 1 : Ly > 

220  DATA  E-6,E-5,E-4, .001, .002, .005, .01, .02, .03 

230  DATA  . I, .2, .3, .4, .3, .6, .7, .3, .9, .95, .93, .99 
240  DATA  .993, .993, .999, .9999, .99999, .999999 

250  READ  /label*1 *> 

260  DATA  l.E-6, l.E-3, l.E-4, .001, .002, .003, .01, .02. .05 

270  DATA  . 1, .2, . 3, .4, .3, .6, .7, .3, .9, .95, .93, .99 
230  DATA  . 995, . 998, . 999, . 9999, . 99999, . 999999 

290  READ  Ycoord' ,*  ■ 

300  ' 

310  N  x  *  1  3 

320  REDIM  Xgr id'  l s  N  ■ 

330  DATA  -7, -6, -5, ~4, -3, “2. -I,  0,1, 2,  3, 4, 5, 6, 7 
340  READ  Xgr i d 1  *  > 

330  1 


360  Nv*  2  7 

370  REDIM  i'gr  i  d ■  1  :  Nt>  1 

330  DATA  l.E-6, l.E-5, l.E-4. , 00 1 , . 002 , . 00 5 , . 0 1 , . 02 . . 05 
390  DATA  . 1 , . 2 . . 3 , . 4 , . 5 , . 6 , . 7 . . 3 , . 9 , . 93 , . 93 , . 99 


400 

DATA  . 

993,  . 993. . 999, . 9999 

4  10 

READ  Y 

gr  id'  *'• 

420 

i 

4  30 

FOP  I* 

1  TO  L" 

44  0 

Y  ■:  o  r  d 

1  I  1  *  F  N  I  n  ■  •  jj  h  i  •  v  c  ■:>  o  r  d 

4  50 

NEXT  I 

4  60 

FOR  I* 

1  TO  N" 

4  70 

•  gr  i  d  ■ 

I  ■  *FN  I  rv-ph  \  •  .  gr  i  d'  I 

4  8  0 

NEXT  I 

4  90 

'.1  =  Ygf 

id-  1  ■ 

500 

•'2*  i gr 

id'll"’ 

E  -3 
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510  GINIT  180.-240. 

520  PLOTTER  IS  505,  "HPGL‘* 

530  PRINTER  IS  505 

540  LIMIT  PLOTTER  505 , 0 . , 1 30 . , 0 . , 240 . 
550  VIEWPORT  20. , 120. , 19. , 132. 

560  '  VIEWPORT  22. , 85. , 59. , 122. 

570  '  VIEWPORT  22. ,35. , 19. ,62. 

580  WINDOW  X 1 , X 2 , V 1 , V 2 

590  '  PRINT  ” V  S  2 “ 

600  FOR  1*1  TO  Nx 

610  MOVE  Xgr  t  d  < I > , Y 1 

620  DRAW  Xgr i d < D , Y2 

630  NEXT  I 

640  FOR  1*1  TO  Ny 

650  MOVE  X 1 , Ygr l d < I  > 

660  DRAW  X2 , Ygr i d< I > 

670  NEXT  I 

680  PENUP 

690  LDIR  0 

700  CSI ZE  2.3,  .5 

710  LORG  5 

720  Y»Y1 -< Y2-Y1 >*. 02 

730  FOR  1*1  TO  Lx 

740  MOVE  Xcoord*'  I ) ,  Y 

750  LABEL  XI ab4 I *■ !• 

760  NEXT  I 

770  CSIZE  3. , . 5 

780  MOVE  .  5  +  •■  X  1  ♦  X 2  ’  ,  V  1  -  .  06 ■  Y  2  -  Y  1  • 

790  LABEL  A* 

300  CSIZE  2. 3, . 5 

310  LORG  3 

320  X*X1 -• X2-X1 >*. 01 

330  FOR  1*1  TO  L" 

340  move  vcoord’  I  • 

350  LABEL  VI ib*l I ■ 

360  NEXT  I 

370  LDIR  PI  2. 

330  C3IZE  3.,. 5 

390  LORG  5 

900  MOVE  :  :i  - .  15-*1  ::2-M  5»*  ,  1  ♦  .'2  • 

910  LABEL  B* 

920  PRINT  " V 3 3 6  " 


*  VERTICAL  PAPER 

!  1  GDU  *  2  mm 

!  TOP  OF  PAPER 
>  BOTTOM  OF  PAPER 


£  4 
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930  Dx*<X2-Xl >'100. 

940  Fl*SQR(Ns'(Ns+l >> 

950  F2*SQR<N*'<Ns*N*-l >> 

960  FOR  1*0  TO  12 

970  Ds*I  '  DEFLECTION  d 

980  Dp-Ds*Fl  !  d£ 

990  FOR  lt»0  TO  100 

1000  T*X  1 +Dx* I t  1  THRESHOLD  T 

1010  Tp*T*F2  '  T/ 

1020  IF  Dp>0.  THEN  1080 

1030  IF  N*  MODULO  2*0  THEN  1060 

1040  Pd*FNPf357<N», Tp> 

1050  GOTO  1120 

1068  Pd*FNPf246<Ns, Tp > 

1078  GOTO  1120 

1088  IF  N*  MODULO  2*8  THEN  1110 

1098  Pd-FNPd357cNs, Tp, Dp> 

1108  GOTO  1120 

1110  Pd*FNPd246< Ni, Tp, Dp ) 

1120  IF  Pd<  *0.  THEN  1160 

1130  IF  Pd>*l .  THEN  1160 

1148  Y*FN I rwph i <  Pd  > 

1158  PLOT  T , Y 

1168  NEXT  It 

1178  PENUP 

1130  NEXT  I 

1190  PAUSE 

1208  PRINTER  IS  CRT 

1210  PLOTTER  505  IS  TERMINATED 

1220  END 

1230  ' 

1240  DEF  FNIrwphi  ..X,«  '  AMS  55,  26.2.2  3 

1258  IF  X* , 5  THEN  RETURN  0. 

1268  P*MIN<  X,  1 .  -:<> 

1270  T*-LOG<P> 

1230  T  *SQR  »T  +  T) 

1290  P»l , *T*< 1 . 4 32733+T- - . 1 39269+T * . 081 38? 

1300  P*T-«  2. 515517  +  T-i  .S02S53-*-T*.  O10323  •  •  F 

1310  IF  .  5  then  p*-p 

1320  RETURN  P 

1330  FNEND 

l  340  ' 

1350  DEP  FNPf 246' DOUBLE  N, PEAL  Tp ■  1  N=2.4.6.... 

1360  DOUBLE  Kj  1  INTEGER 

1370  Pf *. 5-ATN. Tp)  PI 

1330  IF  N*2  THEN  RETURN  Pt 

1390  X*  1  .  <  l  .  *T fjjTp  • 

1 400  S*B  *  l  . 

1410  FOR  )i=l  TO  N  2-2 

1420  B  *B  *:>t  •  l  -  t.V 

1  4  30  S*3-*-B 

1440  NEXT  1} 

1  450  P*'=Pt  -Tp*  » 3  PI 

1460  RETURN  Pt 

1470  FNEND 

14  30  1 


£  -5 
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1490 

DEF  FNPf357< DOUBLE  N,RE«L  Tp>  1 

N*3, 5, 7,  . 

1500 

DOUBLE  Ks  1 

INTEGER 

1510 

X*1 . / < 1 . *T p*T  p > 

1520 

S*flx» 1 • 

1530 

FOR  K s* 1  TO  <N-3)/2 

1540 

flx*0x#X*<Ks-.5)''Ks 

1550 

S*S>0x 

1560 

NEXT  K* 

1570 

Pf  «.  5*<  1 .  -Tp*SQRi'X)  *S  > 

1580 

RETURN  Pf 

1590 

FNEND 

1600 

i 

1610 

DEF  FNPd246<.  DOUBLE  N.RE0L  Tp.Dp.'  ' 

N*2 , 4 , 6 , . 

1620 

Error*l.E-15  1 

ToUrincs 

1630 

Nt*rms*500  1 

Number  of 

1640 

DOUBLE  Ks  ' 

INTEGER 

1650 

X  ■  1 .  •'  •'  1 .  ♦  T  p  *  T  p  > 

1660 

Rk  *30R  C  X  > 

1670 

T  sq*T  p*Rk 

1  630 

D*q*Dp*Dp 

1690 

R*SQR<2.  PI > 

1700 

81-M0-EXP'  5*Dsq  ■  PI 

1710 

8*0 1 +Dp  R 

1720 

B 1  *  .  5*P I -8TN> Tp  • 

1  730 

B-Rk 

1740 

Pd*0  1  *  B  1  ♦0*B 

1750 

FOP  Ks*2  TO  N t  « r Hi i 

1760 

F»FLT(K*-l  ■ 

1770 

T*D.*q*01  F 

1730 

81*3 

1790 

0*  T 

1300 

Rk *Pk *Tsq 

1310 

r».  Pk ♦F  +  Bl  <  k  i 

1320 

B  1  *B 

1330 

B-T 

1340 

P*8  +  B 

1  350 

Pd*  Pd  +  P 

1360 

IF  P  *Err*r-»Pd  THEN  1390 

1370 

NEXT  K  s 

1  330 

PRINT  "500  TERN 3  IN  FNPd246" 

1  390 

IF  N * 2  THEN  PETijPN  Pd 

l  300 

G  1  *  T  i  q  *  E '  P  ■  -  .  5  •  D  i  q  -  •  *  F  N  p  »•,  i  D  p  ■*  T  i  | 

1910 

G*  Tp  Ho  *P  *  Dp •G  1  ■ 

1  920 

F  1  *  Tp  ■*  Dp  P 

1  9  30 

F*Tp*Dp*P 

1  940 

P  d  *  P  d  -  G 

1  950 

IF  N*4  THEN  RETURN  Pd 

1  i60 

FOP  ►»*:  TO  N-i 

l  970 

R  *  ■  t  i  -  1  '  hi 

1930 

T*F 1 *R 

1  990 

F  l  *  F 

2  00O 

F  *  T 

1010 

T  *  *  ■  F  *  G  ♦  R  »  G  1  ■ 

-020 

G  1  *6 

_  0  2  0 

G*  T 

-1  0  4  0 

IF  >■  MODULO  1=1  THEfi  d  =  R  d  -  1 

-  0  j  0 

ME  T  1  : 

_  M  r"  kj 

RETURN  F  j 

,  0  "  0 

FNEND 

1  0  i  0 

1 

E  b 


/  >  v.v.  a’a  /.  /.V/ /.V  ’ •- 
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2090  DEF  FNPd337( DOUBLE  N , REAL  Tp,Dp>  <  N-3,3,7,... 

2100  DOUBLE  K*  '  INTEGER 

2110  X*1 . /( 1 . >Tp*Tp;> 

2120  T*q-Tp*SQR<X> 

2130  D2* . 3*Dp*Dp 

2140  G*T*c|*EXP  < -D2*X  )  *FNPh  l  <Dp*Tso(> 

2130  Pd*FNPh i <  Dp  > -G 

2160  IF  N*3  THEN  RETURN  Pd 

2170  R*SQR  <  2 . ^P I > 

2130  Gl-G 

2190  G»Tp*X*(EXP< -D2  •  P  I  ♦R*Dp*G 1 > 

2200  F 1 *Tp*Dp-R 

2210  F*T  p  +  Dp*R 

2220  FOR  K*«2  TO  N-3 

2230  R* ( K  » -  1  i--Ki 

2240  T*F 1 *R 

2230  Fl-F 

2260  F*T 

2270  T*X*'.  F*l*R*l1  > 

2230  LI *L 

2290  G-T 

2300  IF  Ki  MODULO  2*0  THEN  Pd-Pd-G 

2310  NEXT  K* 

2320  RETURN  Pd 

2330  FNEND 

2340  1 

2330  DEF  FNPh  i  >  'X  •  •  HART,  pig*  140.  *370:?  *3723 

2360  Y*ABS' X '♦.707106731 13634746 

2370  'SELECT  V 

2330  CASE  -  3 . 

2  390  P*  l  6  3  1  .  76026373  37 1 4 70+  i  ♦  ■  4  36 . 26  1  43 3 706 O 926  3  1  ♦  .  ■»  •  36 . 032  "6 22  1  1  94  3  3  :* 
1  0 . 06 4 3339 74 9093423 ♦  i' ♦  .  364  1 ; 39336 76  1  3  1  36  1  4  •  •  ■ 

2  4  00  P*  372  3. 307981 33430672  ♦  .'♦•  7  i  l  3  .  66  32469340493  7*  ,  *•  6733.  2  t  6  -  6  4  l  104?  3 
4  0  3  2 . 2  6  7  0  1  0  8  3  0  0  4  9  7  4  ♦  Y  ♦  P  >  •  • 

24  10  0* 7342 . 4793 1 0 1 9  34  7376*  Y *  1 2 9 6 i . 0 0 4 9 0 14  32  3 0  3 ~ 2 *  -  3 1 ~ . 622  3  3  6  3  043 4 40 

13  3. 077710  730362216* '>'♦'  1  7 . 3  394984  3?  1  393363*  v  •  1  •  ■ 

2420  0*  37  2  3 . 30798 1 33480634 ♦Y* •  1  1  3 1 3 . 1 9208 1  3344033 ♦  -•  13:02. 3 ; 3  -  1 1 4  0  2  0  4 

1  3  349.  34636 1 2  344374* Y*0  ■  •  • 

2  4  30  P  h i * . 3 ♦ E  X  P •  - Y ♦ Y  * ♦ P  0 

2440  C ASE  26 . 6 

2  4  30  P  *  2 . 9  7 :?  3  6  3  6  2  6  3: 9  3  9  9  2  8  9  ■  7 . 4  O  9  7  4  0  6  O  3  9  6  4  "*  4  1  ~  ?  *  ♦  6  .  6  0  2  0  ■* :  3  3  1  3'  3  6  3  0 

3.O1904  *726784  26  746*  <*•  1 . 273  36664472996393»  .  *  .  364  1  3  =•3  3  3  3  4  ~  “330  "  4 
24  60  0*3.  3  6  9  O  7  3  2  O  6  9  3  2  7  3  2  7  7 ♦ i • ■  9 . 6  0  3  9  6  3  3 2  7 1 9 2  7 3  7 3 ~ ♦  *  i ” . o  3 ; 4  4  0  ~  4 “ 4  6  6 0  0 

1  2 . 04 i 93  1  92 78 33  1  2  90 ♦  *  ■  9  .  396 O  3 4 0  l  6 2  3303 4  1  3  ♦ *  2 . 2 6 032 2  3 2 0  ~6  ~  3  2  6  -  ~  * 

2470  P h  1  *  .  3 * E : : P •  -  1' • ■  *  P  Q 

2  4 O  C  A  6  E  ELSE 

2  4  9C  Phi *0 . 

2300  END  SELECT 

2310  IF  0.  THEN  F'h  t  =  1  . -Pi,  1 

-320  RETURN  Phi 

-330  FNEND 
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Table  E-2.  EDF  -  Extreme 


10 

Ns* 1 6  '  NUMBER 

OF 

samples 

20 

Rs*l.  '  RATIO  r 

30 

XI  — 13 

40 

X2*3 

30 

Gun*.  577213664902 

60 

Cl— Gamma 

70 

C2-PI/SQR<6. ) 

80 

F*SQR< .5/PI ) 

90 

fl— 8.  1  LIMITS 

ON 

100 

B*8 .  !  INTEGRAL 

1  10 

Maan  #uiw*0. 

120 

Var  mucw*l,.Ni 

130 

Maan  si gcv*i . 39'Ns 

140 

Var  s  i  gc  v*  1 . 05  1  Ns  + 1 . 5 '• 

150 

Rho*- . 55 

160 

COM  H1,H2 

170 

DIM  A*C301,B*C301 

130 

DIM  XI  aba  1  #<  l  *  30  > ,  V 1  aba  1  * < 1 :  30  > 

190 

DIM  XcoordC 1:30', Ycoord< 1:30) 

200 

DIM  Xgr i  d  < 1 : 30  > , Ygr i d  <  1 : 30  > 

213 

DOUBLE  Lx,Lv,Nx,Ny, I,Ns, It 

220 

1 

230 

AS* " Thr asho 1 d “ 

240 

BS*"Excaadanca  Probabi  1  i  tv" 

250 

1 

260 

Lx*  1  1 

270 

REDIM  XI  ab  a  1  #  '•  1 :  L *  ) ,  Xcoord-  1 :  L  -  • 

230 

DATA  -15, -13, -11, -9, -7, -3, -3,-1, 

1  . 

3,5 

290 

READ  X 1  aba  1  1  ■'  *  > 

300 

DATA  -13, -13, -11, -9, -7, -5, -3.-1, 

1  . 

3.5 

310 

READ  Xc  oor d  < *  > 

320 

1 

3  30 

L  v  *  2  7 

340 

P  E  D I M  Y 1  ab  a  1  $  •  ULv^Ycoord1  t :  L  y  1 

358 

DATA  E-6,E-5,E-4, .001, .002, .005, 

.  0 1  ,  .  0 

2  ,  .  0 

3*0 

DATA  . 6 , . 7 , . 3 , . 9 , . 95 , . 93 , . 99 

370 

DATA  . 995 , . 993 , . 999 , . 9999 , . 99999 

%  • 

Ct  >3  -5 

r*  ’? 

338 

READ  Y  1  aba  IS'*. 

390 

DATA  1 . E-6, 1 . E-5, 1 . E - 4 , . 001 , . 002 

a  • 

0  O  5 , 

.  01  . 

4  00 

DATA  . 6, . 7, . 8, . 9, . 95, . 98, . 99 

4  10 

DATA  .995, .993, .999. .9999, .99999 

a  • 

■Zi  'A  'A  -A 

‘A  ‘A 

420 

READ  Ycoord1 *' 

4  30 

1 

440 

M>  *  1  1 

450 

PEDIM  Xgr id' 1 : N  ■ 

4  <i  0 

DATA  -15, -13. -11. -9, -7, -5. -3,-1, 

l . 

-.  C 
-•  a  J 

470 

READ  ..'grid'  *  1 

4  30 

1 

4  90 

U<’*2  7 

5  0  0 

PEDIM  .grid'  1  :()■•  • 

5  1  0 

DATA  l . E - 6 .  1  . E - 5 .  1  . E - 4 .  .001,  .002 

,  , 

0  m  5  . 

.  0  1  , 

520 

Z  *  M 

DATA  ,  j'ic.  ,  .  ,  .  i-i-s,  ,  iiSi,  .  S'ji'V; 

M  ‘A  -Z*  Zt 

U  3 4 

5  4  0 

READ  .  gr  id  *' 
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FOR  1*1  TO  Ly 

YcoordC  I  )*FNIn>jphi  <Ycoord<  I  >  > 

NEXT  I 

FOR  I-i  TO  Ny 

Ygr 1 d< I  >*FNIr»ophi  <Ygrid<I>> 

NEXT  I 
Yl*Ygr i d<  1  > 

Y2*Ygr i d(Ny ) 

GINIT  180.^240. 

PLOTTER  IS  305, “ HPGL “ 

PRINTER  IS  503 

LIMIT  PLOTTER  505, 0. , 130. , 0. , 240. 
VIEWPORT  20. , 120. , 13. , 132. 
VIEWPORT  22. ,35. ,53. , 122. 

VIEWPORT  22. ,35. , 13. ,32. 

WINDOW  Xl,X2,Yi,Y2 
PRINT  "VS2" 

FOR  1*1  TO  Nx 
MOVE  Xgr i d< I > , Y 1 
DRAW  Xgr l d< I > , Y2 
NEXT  I 

FOR  1*1  TO  Ny 
MOVE  XI ,  Ygr  id1;  I  • 

DRAW  X2  ,  Ygr  i  d  <.  I  ' 

NEXT  I 
PENUP 

CS  1 2E  2. 3, . 5 
LORG  5 

Y  *  Y  1  - 1  r'2-Y  1  >*.02 
FOR  1*1  TO  Lx 
MOVE  Xc oor d < I > , V 
LABEL  XI  ib«  1  $>■  I  ' 

NEXT  I 
CSIZE  3. , . 5 

MOVE  .  5*  1  X  1  *::2  • ,  V  1  -.  03*  -  r2-  Y  1  • 
LABEL  A* 

CSIZE  2.  3,  .5 
LORG  3 

x*x l  - •  :;2-: :  l  •*.  Oi 

FOP  1*1  TO  L" 

MOVE  X  ,  V.;  oor  d  ■  I  1 
LABEL  1  ib*  1  *  <  I  • 

NEXT  I 
LDIP  PI  2. 

CSIZE  3... 5 
LORG  5 

HO-'E  1  -  .  1  5  *  .  2  -  1  •  .  .  5-  .  1  ♦  ■  2 

LABEL  Bf 

PRINT  33" 

D  = 1  2-1  lOO. 

D  I  M  P  0:1  0  0 


>  vertical  paper 

1  1  G DU  *  2  mm 

1  TOP  OF  PAPER 
1  BOTTOM  OF  PAPER 


► 


► 


I 


\ 

i 

i 
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1070  FOR  D**0.  TO  7.5  STEP  .75  !  DEFLECTION  d 

1080  FOR  lt*0  TO  100 

1090  J-It 

1100  Th«X2~Dx*It  '  THRESHOLD 

1110  Mut»Mean_muco+Th*f1ean_si  gcu 

1120  T  «mp*Var_muc  y  +  Th*Th*Var_s i  gc  v 

1130  T«mps*T*mp+2 .  #Th*SQR<  Var_muc  y*Var_ s  i  gc  u )  *Rho 

1140  Si gmat*SQR<T«mp) 

1150  Hl«Cl+C2*<Mut-Ds)/'Rs 

1160  H2»C2*Si gmat/Rs 

1170  Sa«FNS<fl> 

1180  Sb*FNS<  B  > 

1190  PRINT  “FNS(ft>  =  “ ;Sa;"  FNS(B>  »  “ ; Sb 

1200  DOUBLE  N,K,J 

1210  N»2 

1220  H* < B-fl )  * .  5 

1230  $*< Sa+Sb>*. 5 

1240  V-1.E10 

1250  T»0. 

1260  FOR  K*  i  TO  N-l  STEP  2 

1270  T»T+FNS<fi+H*K) 

1280  NEXT  K 

1290  S*S+T 

1300  Vo*V 

1310  V«<S+T>*H*2. /3. 

1320  IF  ABS<V-VoX  1 .  E-9  THEN  1380 

1330  N-N+N 

1340  H»H*.5 

1350  GOTO  1250 

1360  Pd*F*V 

1370  IF  Pd >.9999995  THEN  1400 

1380  P<  It  )*FNIrwphi  <Pd> 

1390  NEXT  It 

1400  FOR  1 1  *0  TO  J-l 

1410  PLOT  X  2 - D  x  *  1 1 , P ■  It  > 

1420  NEXT  It 

1430  PENUP 

1440  NEXT  Ds 

1450  PAUSE 

1460  PRINTER  IS  CRT 

1470  PLOTTER  505  IS  TERMINATED 

1430  END 

1490  ! 

1500  DEF  FNIrvphi  •. X  '  '  RMS  55,  26.2.23 

1510  IF  X=.  5  THEN  RETURN  0. 

1520  P  *  M  I  N »  X ,  1 . - X  1 

1530  T»-LOGvP^ 

1540  T»SQRvT+T> 

1553  P=1.*T»'  1.432733  +  T*-  .  1 39269  +  T-*  .  OO  1  308  •  • 

1560  P*T-'  2. 515517  +  T  —  .  802353  + T-r .  010323  •  ■  P 

1570  IF  X  .5  THEN  P*-P 

1530  RETURN  P 

1590  FNEND 

1600  ' 

16  10  DEF  FNS' X  1 

1620  COM  H1.H2 

16  30  T*.  5-x+x+exf-  hi  +  h2*;;  • 

1 6  4  i_t  IF  T  ICO.  THEN  RETURN  0 . 

1650  RETURN  E:.R'-T' 

1660  FNEND 
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Table  E-3.  Simulation  -  Extreme 


10 

Ns*  1 6 

| 

NUMBER  OF  SAMPLES 

20 

DATA  0, .75, 1.5,2.25, 

3. ,3. 75,4 

. 5, 5. 25, 6. , 6. 75, 7. 5  !  DEFLECTION 

30 

Rs*l . 

i 

SCALING  r 

40 

Nt*2~23 

! 

NUMBER  OF  TRIALS 

30 

A*=" E- 1 6-75 to75- 1-22 

■I  i 

FILE  NAME 

60 

Nb*1000 

i 

NUMBER  OF  BINS 

70 

Mean*-. 577215664902 

I 

FOR  Vt,  LOG  OF  EXPONENTIAL 

30 

Si gma*PI 'SQR<6.  > 

j 

RANDOM  VARIABLE 

30 

Zmi n*-15. 

! 

MINIMUM  Z  VALUE 

100 

Zmax*5. 

f 

MAXIMUM  Z  VALUE 

110 

DOUBLE  Ns,Nt , Nb, Kt , K 

s ,  J  • 

INTEGERS 

120 

DIM  P0C 1000 >, PIC  1000 

> , P2< 1000 

>  ,  P3  <  1 000  >  ,  P4 1  000  >  ,  P5  1  000  ' 

130 

DIM  PS  <  1  000  >  ,  P7  (.  1  000 

),P8<1000 

> ,  P9C 1 000  > ,  P 1 0 '  1 000  > ,  Edf 1:  1  1  000  > 

140 

READ  D0,D1,D2,D3,D4, 

D  5 ,  D  6  ,  D  7 , 

D3, D9, D10 

150 

A0=-Me  an /  S  i  gma 

160 

Al  =  l . /Si gma 

170 

F  =  1 .  /  N  s 

130 

G»1 . /SGR  <  Ns- 1  > 

190 

Dz*<Zmax-Zmi n>/Nb 

200 

FOR  Kt  *  1  TO  Nt 

1 

SIMULATION 

210 

31*32=0. 

220 

FOR  Ks*l  TO  Ns 

230 

V  t  *  L  0  G  <  -  L  0  G  <  R  N  D  ■'  > 

i 

UN-NORMALIZED  RANDOM  VARIABLE. 

240 

V=A1 *Vt +A0 

i 

ZERO  MEAN,  UNIT  VARIANCE  RV 

250 

31*31 +  V 

260 

S2=S2+V*V 

270 

NEXT  Ks 

230 

M  c  *  3 1  *  F 

i 

SAMPLE  MEAN 

290 

3  c  *  3  Q  R  1 3  2  -  N  s  *  M  c  *  M  ■:  >  * 

G  ! 

SAMPLE  STANDARD  DEVIATION 

300 

Vt  =L0G  <  -LOG  (  PND  .>  > 

310 

V*A 1  * Vt  +A0 

320 

C=Rs*V-Mc 

330 

Z*1  D0  +  C  3c 

340 

J  =  I N  T  •:  1  Z  —  Z  n  i  i  n  > D  ’  ) 

350 

IF  J  0  THEN  J  *  0 

360 

IF  J  Nb  THEN  J=Nb 

370 

P  0  ‘  J  1  *  P  0  <  J  1  ♦  1  . 

330 

Z  *  ■  D 1  +  C  1  3  •: 

390 

J  =  I  N  T  1  k  Z  -  Z  (nt  n  >  D  z  > 

4  0  0 

IF  J  0  THEN  J*0 

4  1  0 

IF  J  Nb  THEN  J=Nb 

420 

PI'  J  '  *P  1  ■  J  '>  ♦  1  . 

4  30 

Z*  > D2+C  >  3c 

44  0 

J  *  I  N  T  ■  i.  Z  -  Z  fti  i  r,  ■  D  z  ■ 

450 

IF  J  O  THEN  J  *  0 

4  60 

IF  .1  Mb  THEN  J  =  H  fc< 

4  70 

P2 ■  J  ■  =  F 2 '  J  '*1  . 

4  3  0 

Z *  ■  D  3 ♦  C  ■  3c 

4  9  0 

J  =  I  NT  >  •  Z -Zru  i  n  ■  D z  • 

5  O  0 

IF  I  O  THEN  1=0 

5  1  0 

IF  r  Mb  THEM  J  =  Mb 

520 

P  3  •  J  •  *  P  3  ■  J  ■  ♦  1  . 
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330 

2*<D4+C)/Sc 

540 

J»INT<(2-2min>^Dz) 

530 

IF  J< 0  THEN  J»0 

560 

IF  J>Nb  THEN  J°Nb 

570 

P4  <  J  >  *P4 <  J )  + 1 . 

580 

2*<D3+C)/Sc 

590 

J«INT<  <2-2mi n)/Dz) 

600 

IF  J<0  THEN  J=0 

610 

IF  J>Nb  THEN  J=Nb 

620 

ps< J)*P3< J>+i . 

630 

2*<D6+C>/S c 

640 

J»INT<  <2-2mi n)/Dz> 

650 

IF  J<0  THEN  J=0 

660 

IF  J >Nb  THEN  J  =  Nb 

670 

P6  <  J ) *P6<  J )  + 1 . 

630 

2*<D7+C)',Sc 

690 

J»INT(<Z-Zmin>/Dz> 

700 

IF  J<0  THEN  J»0 

710 

IF  J >Nb  THEN  J  =  Nb 

720 

P7< J)*P7< J)+l . 

730 

2*<D8+C)^  Sc 

740 

'  J*INT<  <Z-Zmi  rO'  Dz) 

750 

IF  J<0  THEN  J=0 

760 

IF  J >Nb  THEN  J=Nb 

770 

P8< J)«P8< J>+1. 

730 

2*<D9+C)/Sc 

790 

J*INT < (2-2mi n>  Dz) 

300 

IF  J<0  THEN  J=0 

310 

IF  J  >Nb  THEN  J=Nb 

320 

P9<  J )*P9<  J> ♦ 1 . 

330 

2a<D10+C>/Sc 

340 

J= I NT <  <  2-2m 1 n > • Dz ) 

350 

IF  J<0  THEN  J=0 

360 

IF  J  >Nb  THEN  J«Hb 

370 

P  1  0  >  J  >  *  P  1  0  <  J  >  ♦  1  . 

330 

NEXT  Kt 

390 

MAT  P0*P0  -  Nt  > 

900 

MAT  P1*P1  •  Nt  > 

910 

MAT  P2=P2  '  < Nt  > 

920 

MAT  P3*P3  -  <•  Nt  > 

930 

MAT  P4*P4/<Nt) 

940 

MAT  P 5 » P 5  <  N  t  > 

930 

MAT  P 6 * P 6 '  <  NO 

960 

MAT  P?=P7  •:  Nt  > 

970 

MAT  P 3 a P 3  '  N X  '> 

980 

MAT  P9  =  P9  '  Nt  '■ 

990 

MAT  P10=P10  > Nt  • 

2 
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8  Jun  198?  14: 18:02 

.  1010 

FOR  J*Nb  TO  0  STEP  -1 

1020 

S0«S0+P0< J> 

1030 

IF  S0-0.  THEN  1060 

1040 

IF  S0>»1.  THEN  1070 

1050 

P0< J)*FNInvphi <S0>  !  P0< J>*Prob<Z0>*Zmin+J*Dz> 

1060 

NEXT  J 

1070 

P0< J>*0. 

1080 

FOR  J«Nb  TO  0  STEP  -1 

1090 

S1-S1+P1 < J) 

1100 

IF  S1*0.  THEN  1130 

1110 

IF  S1>»1.  THEN  1140 

1120 

PI  < J>*FNInvphi <S1 ) 

1130 

NEXT  J 

1140 

Pl< J>»0. 

1150 

FOR  J*Nb  TO  0  STEP  -1 

1160 

S2*S2+P2  <  J  > 

1170 

IF  S2*0.  THEN  1200 

1130 

IF  S2>*1.  THEN  1210 

1190 

P2< J)*FNInvphi <32> 

1200 

NEXT  J 

1210 

P2  <  J  >  *0 . 

1220 

FOR  J  =  Nb  TO  0  STEP  -1 

1230 

S3*S3+P3  <  J  ) 

1240 

IF  33*0.  THEN  1270 

1250 

IF  S3>*1.  THEN  1230 

1260 

P3<  J>*FNIriMphi  <S3> 

1270 

NEXT  J 

1280 

P3< J)»0. 

1290 

FOR  J  =  Nb  TO  0  STEP  -1 

1300 

S4*S4+P4  <  J ) 

1310 

IF  54=0 .  THEN  1340 

1320 

IF  S4  >■ l .  THEN  1350 

1330 

P4  < J > *FN I nvph  t <  3  4  > 

1340 

NEXT  J 

1350 

P4( J)«0. 

1360 

FOR  J* Nb  TO  O  STEP  -1 

1370 

35=35 +P5 1  J  > 

1  380 

IF  35*0.  THEN  1410 

1  390 

IF  3  5  >  *  1  .  THEN  1420 

1400 

P5<.  J  >  *FNI  nvph  l  i  35. > 

1410 

NEXT  J 

1420 

P5  <  J  >  *0 . 

14  30 

FOR  J  =  Nb  TO  0  STEP  -1 

1440 

3  6  *  8  6  +  P  6 '  J  > 

1450 

IF  86*0.  THEN  1480 

1  460 

IF  3  6  >  =  1 .  THEN  1490 

14  70 

P  6  <.  J  >  *  F  N  I  Ci  • 1  p  h  l  •  8  6  1 

14  80 

NEXT  J 

1490 

P6  V  J  '  *  0 , 

1500 

FOR  J=Hb  TO  0  STEP  -1 

1510 

37  =  87  +  P  7 ■  ]  ■ 

1520 

IF  37=0.  THEN  1550 

15  30 

IF  87  =1.  THEN  1560 

154  0 

P7  '  T  •=FNIrvph  i  '  37  ■ 

1  550 

NEXT  J 

1  560 

P 7  '  J  -30. 
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1370  FOR  J»Nb  TO  0  STEP  -1 

1380  S8-S8+P8<J> 

1390  IF  S8-0.  THEN  1620 

1600  IF  S8>«1.  THEN  1630 

1610  P8< J>*FNInvphi <S8> 

1620  NEXT  J 

1630  P8<J>»0. 

1640  FOR  J*Nb  TO  0  STEP  -1 

1650  S9»S9+P9<J> 

1660  IF  S9»0.  THEN  1690 

1670  IF  S9  >* 1 «  THEN  1700 

1680  P9< J>*FNIrwphi <S9) 

1690  NEXT  J 

1700  P9<J)-0. 

1710  FOR  J»Nb  TO  0  STEP  -1 

1720  S 1 0*S 1 0+P 1 0  <  J  > 

1730  IF  S10»0.  THEN  1760 

1740  IF  S10>»1.  THEN  1770 

1 750  P 1 0  ( :  J ) *FN I rwph i <  S 1 8  > 

1760  NEXT  J 

1770  P10<J)»0. 

1730  FOR  J-l  TO  Nb 

1790  Edf < J)*P0< 

1800  Edf  <  J+Nb)*Pl  J> 

1318  Edf < J+Nb*2>*P2< J> 

1320  Edf ( J+Nb*3 > *P3<  J > 

1330  Edf  <  J+Nb*4  >*P4<  J> 

1340  Edf < J+Nb*3>*P5< J) 

1350  Edf ( J+Nb*6  >  *P6< J  > 

1360  Edf<J+Nb*7>«P7<J> 

1370  Edf < J+Nb*8>*P8< J> 

1330  Edf < J+Nb*9>*P9< J > 

1390  Edf  <  J+Nb*  1 0  >  *P  1 0  <:  J  > 

1900  NEXT  J 

1910  CREATE  DATA  A*, 396 

1920  ASSIGN  #1  TO  A* 

1930  PRINT  *ljEdf<*> 

1940  ASSIGN  #1  TO  * 

1950  PAUSE 

1960  END 

1970  ' 

1930  DEF  FN I iwph  i  X  >  1  AMS  55,  26. 2.  23 

1990  IF  X * . 5  THEN  RETURN  0. 

2000  P=M I N  <  X , 1 . -X  > 

2010  T  =  - L 0 G 1 P  > 

2020  T =SQR •:  T  +  T  > 

2030  P=  l .  -t-T*  ■  l .  432733+T*-  .  1  39269  +  T*  .  OO  1  303  •  - 
2040  P=T--  2. 515517+T--  .302S53  +  T*. 010323  ■  •  P 

2050  IF  X  . 5  THEN  P  =  -P 
2060  RETURN  P 

2070  FNEND 
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Table  E-4.  Plot  -  Simulation 


A**"E-l6-73to73-i-23" 

Zmi n*-13. 

Zmax*3. 

Xl»-13. 

X2-3. 

Nb*1000 

ASSIGN  #1  TO  A* 

READ  #1 ; Edf  <*> 

ASSIGN  #1  TO  * 

DIM  Edf < 1 s 11000) 

DIM  A#C30] , B*C303 

DIM  X 1  abe 1  *  < 1 : 30 ) , Y 1 abe 1  *  <  1 : 33 ) 

DIM  XcoordC 1 s  30) , Ycoord< 1:30) 

DIM  Xgri d< 1 : 30) , Ygr i d< 1 : 30) 

DOUBLE  Nb, Lx, Ly, Nx, Ny, I , K  !  INTEGERS 

A**“Thresho1 d“ 

B**“ Exceedance  Probabi  1  i  t  v" 

Lx*  1 1 

REDIM  X 1  abe  1 1  <1 :  Lx  ) ,  Xc  oord  <  1 :  Lx  ;• 

DATA  -13, -13, -11, -9, -7, -5, -3,-1, 1,3,5 
READ  X label *<*) 

DATA  -13, -13, -11, -9, -7, -5, -3, -1,1, 3, 5 
READ  XcoordC*) 

j 

Ly-27 

REDIM  Yl4beI#<l:Ly>,Ycoord<l:Ly) 

DATA  E-S, E-5, E-4, . 001 , , 002, . 005, . 01 , . 02, . 05 
DATA  .  1 ,  .  2 ,  .  3 ,  .  4 ,  .  5 ,  .  S ,  .  •  8 ,  .  9 ,  .  95 ,  .  93 ,  .  99 

DATA  .  995  ,  .  993 , . 999 , . 9999 , . 99999 , . 999999 
READ  Y 1 ab e 1  * <  * ) 

DATA  1 . E-6, 1 . E-5 , 1 . E-4, . 001 , . 002, . 005, .01, .02, . 03 
DATA  .1,. 2, .3, .4, .5, .6, .7, .3, .9, .95,. 93 , . 99 
DATA  .995, .993, .999, .9999, .99999, .999999 
READ  Ycoord<*> 

! 

Nx*  1*1 

REDIM  XgridCl:Nc> 

DATA  -15, -13, -1 1 , -9, -7, -5, -3, -1 . 1 , 3, 5 
READ  Xgr id' *■ 

i 

Ny*27 

REDIM  Ygr  i  d-  1:  M"  ■ 

DATA  1  . E-8, 1 .E-5,  1 . E-4, . 001 , . 002. . 005. .01,-02.. 0‘ 
DATA  .  1 ,  . 2 . . 3 , . 4 . . 5 . . 3 . . 7 . . 3 , . 9 . . 95 ,  . 98 ,  . 99 
DATA  . 995. . 993. . 999, . 9999, . 99999. . 999999 
READ  Ygr i d 1 « • 
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500 

FOR  1-1  TO  Ly 

510 

Yc  oord  <  I  >=FNIrwph 

520 

NEXT  I 

530 

FOR  I»1  TO  Ny 

540 

Ygr  i  d< I >*FNInvphi 

550 

NEXT  I 

560 

Yl=Ygri d<  1 ) 

570 

Y2*Ygr i d<  Ny ) 

580 

GINIT  180./240. 

590 

PLOTTER  IS  505, ”H 

600 

PRINTER  IS  505 

610 

LIMIT  PLOTTER  505 

620 

VIEWPORT  20. , 120. 

630 

!  VIEWPORT  22., 35., 

640 

!  VIEWPORT  22. ,35. , 

650 

WINDOW  X 1 , X2 , Y 1 , Y 

660 

!  PRINT  " VS2“ 

670 

FOR  I»1  TO  Nx 

630 

MOVE  Xgr i d< I > , Y1 

690 

DRAW  Xgr i d< I > , Y2 

700 

NEXT  I 

710 

FOR  1  =  1  TO  Ny 

720 

MOVE  XI , Ygr i d<  I  > 

730 

DRAW  X2 , Ygr i d<  I  > 

740 

NEXT  I 

750 

PENUP 

760 

LDIR  0 

770 

CSI2E  2.3, .5 

730 

LORG  5 

790 

Y  ■  Y 1  -  < Y  2  -  Y  1  >*.02 

800 

FOR  1=1  TO  Lx 

310 

MOVE  Xc  oor d  <  I  >  ,  Y 

320 

LABEL  XI abel *< I  1 

330 

NEXT  I 

840 

CSI2E  3. , .5 

350 

MOVE  .  5*<:X1+X2>  .  V 

360 

LABEL  A* 

370 

CSIZE  2. 3, . 5 

330 

-ORG  3 

390 

::=xi ->x2-xi  ^*.01 

900 

FOR  1=1  TO  Ly 

910 

MOVE  X, Ycoordv  I  1 

920 

LABEL  Y1  atpel  I  • 

930 

NEXT  I 

940 

LDIR  PI  2. 

950 

CSIZE  3. . . 5 

960 

LORG  5 

970 

move  xi-. 15-'  : : 2 - : 

930 

LA3EL  B* 

990 

FPINT  "VS  36" 

i  VERTICAL  PAPER 


i  1  GDU  =  2  mm 

i  TOP  OF  PAPER 
!  BOTTOM  OF  PAPER 


1  - .  06- 


:2-::i  •  .  .5  —  .-»♦ 
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Dz*<2max-2mir)>/Nb 
FOR  1*1  TO  Nb 
T*Edf <1 ) 

IF  T-0.  THEN  1830 
PLOT  Zm i n  +  Dz* I , T 
NEXT  I 
PENUP 

FOR  1*1  TO  Nb 
T*Edf < I+Nb; 

IF  T*0.  THEN  1110 

PLOT  Zmin+Dz*I,T 

NEXT  I 

PENUP 

K«Nb*2 

FOR  1*1  TO  Nb 

T *Edf  <  I  -HO 

IF  T*0 .  THEN  1130 

PLOT  2mi n+Dz* I , T 

NEXT  I 

PENUP 

K*Nb*3 

FOR  1*1  TO  Nb 
T*Edf < I+K> 

IF  T *0 .  THEN  1250 

PLOT  2m) n+Ds*I » T 

NEXT  I 

PENUP 

K*Nb*4 

FOR  1*1  TO  Nb 

T*Edf < I -HO 

IF  T *0 .  THEN  1320 

PLOT  2m i n+Dz* I , T 

NEXT  I 

PENUP 

K*Nb*5 

FOR  1*1  TO  Nb 
T*£df <■  I +  K  > 

IF  T*0.  THEN  1390 

PLOT  Zm l n+Dz* I , T 

NEXT  I 

PENUP 

K*Nb*6 

FOR  1*1  TO  Nb 
T  *  E  d  t'  <  I  ♦  K  > 

IF  T  =  0  .  THEN  1460 
PLOT  2  rn  '  n  +  D  z  ■*  I .  t 
NEXT  I 
PENUP 
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1480  K*Nb*7 

1490  FOR  1*1  TO  Nb 

1500  T  *Edf  <1+10 

1510  IF  T*0.  THEN  1530 

1520  PLOT  Zmin+Dz*I,T 

1530  NEXT  I 

1540  PENUP 

1550  K*Nb*8 

1560  FOR  1*1  TO  Nb 

1570  T*Edf<I+K; 

1580  IF  T *0 .  THEN  1600 

1590  PLOT  Ztn  i  n+Dz*I ,  T 

1600  NEXT  I 

1610  PENUP 

1620  K*Nb*9 

1630  FOR  1*1  TO  Nb 

1640  T*Edf < I+K> 

1650  IF  T =0 .  THEN  1670 

1660  PLOT  Zmin+Dz*I,T 

1670  NEXT  I 

1680  PENUP 

1690  K*Nb* 1 0 

1700  FOR  1*1  TO  Nb 

1710  T  *Edf  <  I -HO 

1720  IF  T  =  0 .  THEN  1740 

1730  PLOT  Zm i n  +  Dz* I , T 

1740  NEXT  I 

1750  PENUP 

1760  PAUSE 

1770  PRINTER  IS  CRT 

1780  PLOTTER  505  IS  TERMINATED 

1790  END 

1800  ! 

1  o  I  0  D  E  F  F  N  I  r'i 1 1  p  h  i  1  *! .  '  1  AMS  55 

1820  IF  X * . 5  THEN  RETURN  O. 

13  30  P=M I N 1  1 . -X  > 

1340  T *-LQG< P  > 

1350  T *SOR  <  T  +  T  > 

1360  P*  1  .  +  T-*  1  1 . 4  32783+T  *  •  .  1 39269 +T 

1370  P  =  T-<  2. 515517  +  T-*'  .  802S53  +  T-*.  0 

1380  IF  :<•  .  5  THEN  P*-P 

1390  RETURN  P 

1900  FNEND 
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